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1  Summary 


Development  of  the  equations  to  predict  the  species  (H  + ,  CH  3  OH,  H  2  O)  molar  fluxes 

and  ionic  current  density  through  the  PEM  of  a  DMFC  is  described.  In  the  Transport  Set 

1  section  of  this  report,  the  equations  are  based  on  the  extension  of  moderately  dilute 
solution  theory.  These  equations  involve  the  Fick’s  law  type  diffusivities  that  are 
strongly  dependent  on  the  species  concentrations.  In  the  Transport  Sets  II  and  III 
sections  of  this  report,  the  developments  are  valid  for  dilute  to  very  strong  solutions.  The 
transport  set  II  is  the  application  of  the  generalized  Maxwell-Stefan  equations  for  the 
species  transport  through  the  PEM.  The  transport  through  set  III  is  the  result  of  the 
application  of  the  Lars  Onsager’s  irreversible  thermodynamic  approach  to  transport 
processes. 

2  Introduction 

Modeling  activity  in  conjunction  with  the  experimental  work  is  deemed  to  economically 
aid  in  the  development  of  an  optimal  design  of  a  direct  methanol  fuel  cell  (DMFC)  for 
aerospace  and  ground  transportation  vehicle  applications.  To  this  end,  research/ 
development  activity  by  the  principal  investigator  (P.I.),  Sarwan  S.  Sandhu  is  given 
below. 

The  principal  investigator  (P.I.)  developed  a  theoretical  scheme  (see  Volume  1  of  this 
report)  to  predict  the  reversible  cell  voltage  in  the  absence  of  electrode  poisoning  and 
methanol  crossover  through  a  polymer  electrolyte  membrane  of  a  DMFC,  for  example, 
Nation  perfluorosulfonic  acid  polymer  and  polybenzoxazole  or  polybenzimidazole 
based  electrolyte  membranes.  The  scheme  based  on  the  application  of  thermodynamic 
fundamentals  of  the  phase  and  electrochemical  reaction  equilibria  predicts  the  reversible 
voltage  at  any  fuel  cell  temperature.  The  effect  of  nonideal  behavior  of  fluid  phases  in 
the  anodic  and  cathodic  fluid  compartments  of  a  fuel  cell  on  the  reversible  voltage  is 
accounted  for.  The  simulation  is  capable  of  predicting  the  reversible  cell  voltage  as 


1 


function  of  reactant  feed  composition  and  different  total  pressures  in  the  anode  and 
cathode  plenums.  The  developed  reversible  voltage  equation  reduces  to  the  Nernst  type 
equation  for  the  ideal  fluid  phase  behavior  of  the  fuel  and  oxidant  feeds.  The  predicted 
reversible  cell  voltage  can  be  used  as  an  ideal  standard  to  which  the  actual  voltage  of  a 
direct  methanol  fuel  cell  operating  at  a  given  set  of  reactant  feed  composition, 
temperature,  and  pressure  conditions  can  be  compared  to  evaluate  the  fuel  cell 
performance. 

A  senior  year  undergraduate  student,  Mr.  R.  Owen  Crowther,  typed  the  entire  theoretical 
development  by  the  P.I.  for  the  reversible  cell  voltage  prediction.  The  P.I.  guided  Owen 
to  develop  a  computer  code  to  simulate  the  entire  set  of  mathematical  equations  to 
generate  the  reversible  cell  voltage  data  as  function  of  temperature,  fuel  feed 
composition,  and  anode  and  cathode  side  total  pressures. 

The  actual  theoretical  development,  the  computer  code,  and  the  code  generated  numerical 
data  in  the  form  of  tables  and  plots  are  presented  in  Report  #1  entitled:  “Fuel  Cell 
Project:  Direct  Methanol  Fuel  Cell:  Theoretical  Formulation  of  Reversible,  Open-Circuit 
Voltage  for  a  Direct  Methanol  Fuel  Cell”  by  Sarwan  S.  Sandhu,  Department  of  Chemical 
&  Materials  Engineering,  The  University  of  Dayton,  Dayton,  Ohio  45469-0246. 

Presented  in  this  report  volume,  the  P.I.  has  developed  the  mathematical  equations  to 
predict  the  ionic  current  density  and  the  species  (H+,CH3OH,  FEO)  molar  fluxes  for  their 
transport  through  a  solid  polymer  electrolyte  or  polymer-ceramic  material  composite 
membrane  of  a  DMFC  at  the  steady  state,  isothermal  conditions  by  the  application  of 
fundamentals  of  transport  phenomena.  The  developed  transport  equations  are  presented 
in  three  sets.  In  the  transport  set  I,  the  equations  are  based  on  the  extension  of 
moderately  dilute  solution  theory.  These  equations  involve  the  Fick’s  law  type  mass 
diffusivities  that  are  known  to  be  strongly  dependent  on  the  species  concentrations.  The 
transport  equation  sets  II  and  III  are  valid  for  very  strong  solutions.  The  transport 
equation  set  II  is  based  on  the  application  of  the  generalized  Maxwell-Stefan  equations. 
The  transport  equation  set  III  is  based  on  the  application  of  the  Lars  Onsager’s 
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irreversible  thermodynamic  approach  to  transport  processes.  The  developed  equations 
account  for  the  effect  of  the  voltage,  pressure  and  the  species  concentration  gradients  on 
their  transport  fluxes.  These  equations,  describing  species  transport  fluxes  through  a 
solid  polymer  electrolyte  membrane  are  to  be  coupled  with  the  equations,  yet  to  be 
developed,  describing  mass  transfer  through  the  electrode  porous  backing  layers  and 
species  mass  transport  with  electrochemical/chemical  kinetics  in  the  porous  electrode 
layers.  The  coupled  equations,  describing  the  various  phenomena  occurring  in  a  DMFC 
in  operation,  would  be  of  immense  significance  in  the  evaluation  of  perfonnance  and 
design/development  of  a  DMFC.  The  usefulness  of  the  developed  transport  equations 
may  also  be  seen  in  that  they  can  be  employed  to  experimentally  determine  the  transport 
properties  of  a  membrane  such  as  permeability  of  the  fluid  mixture  and  effective  mass 
diffusivities  of  the  species  by  designing  appropriate  fuel  cell  experiments. 

The  developed  transport  equations  by  the  P.I.  were  typed  by  Mr.  Sarath  Krishnan,  a 
graduate  student  in  the  chemical  engineering  program  at  the  University  of  Dayton.  Very 
recently,  Mr.  R.  Owen  Crowther  (who  is  now  a  graduate  student  in  our  chemical 
engineering  program)  has  been  guided  by  the  P.I.  to  simulate  the  first  set  of  transport 
equations  in  the  form  of  a  computer  code  to  generate  numerical  data  on  current  density 
and  species  (H+,  CH3OH,  H2O)  molar  fluxes  through  a  polymer  electrolyte  membrane  of 
a  DMFC.  He  is  expected  to  start  developing  the  computer  code  in  the  near  future. 

The  actual  theoretical  development  of  the  transport  equations  is  presented  in  Report  #2 
entitled:  “Fuel  Cell  Project:  Direct  Methanol  Fuel  Cell:  Theoretical  Formulation  of 
Transport  Fluxes  of  species  (H+,  CH3OH,  H2O)  through  a  Solid  Polymer  Electrolyte 
Membrane  (PEM)  of  a  Direct  Methanol  Fuel  Cell,”  by  Sarwan  S.  Sandhu,  Department  of 
Chemical  &  Materials  Engineering,  The  University  of  Dayton,  Dayton,  Ohio  45469- 
0246. 

The  P.I.  developed  a  set  of  equations  to  predict  the  reversible  as  well  as  irreversible  or 
actual  power  production  from  a  continuously  fed  DMFC  operating  at  steady-state 
conditions.  This  development  is  based  on  the  fundamentals  of  classical  thennodynamics 
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and  the  species  and  overall  material  balances.  The  scheme  requires  the  information  on 
the  reactant  feed  stream  composition,  temperature,  and  pressure  conditions  at  the  inlets  of 
a  DMFC,  fractional  conversion  of  methanol  via  electrochemical  process  for  the 
production  of  electric  power,  its  chemical  oxidation  at  the  catalytic  surface  of  the 
cathode,  and  the  amounts  of  methanol,  water,  and  hydrogen  ions  transported  through  the 
solid  polymer  electrolyte  membrane  from  the  anode  side  to  the  cathode  side  of  a  DMFC. 
The  development  accounts  for  the  nonideal  behavior  of  the  fluid  phases.  This 
development  is  not  in  the  typed  form  at  the  time  of  reporting.  Mr.  Sarath  Krishnan,  a 
graduate  student  in  the  chemical  engineering  program  at  the  University  of  Dayton  has 
gladly  accepted  to  understand  the  entire  scheme  of  equations  and  to  present  them  in  the 
form  of  a  typed  report.  Also,  he  would  develop  a  computer  code  under  my  direction  to 
generate  numerical  data  on  the  reversible  and  actual  electric  power  production  as  function 
of  input  parameter  values. 

The  P.I.  studied  and  carried  out  partial  analysis  of  the  open  literature  on  the  research/ 
development  of  direct  methanol  fuel  cells  depending  on  availability  of  time  for  this 
equally  important  activity.  Based  on  the  information  acquired  from  the  literature  on  the 
direct  methanol  fuel  cells  and  his  knowledge  of  the  basic  and  engineering  sciences 
relevant  to  the  research/development  of  a  DMFC,  the  P.I.  developed  an  outline  of  a 
comprehensive  research/development  program  to  develop  an  efficient  DFMC  energy 
conversion  system.  The  outline  has  been  entitled,  “Direct  Methanol  Fuel  Cell  (DMFC) 
Research/Development  -  Theoretical  and  Experimental  Integrated  Approach”.  A  copy  of 
this  outline  was  provided  to  Dr.  J.P.  Fellner  of  the  Power  Division  of  the  WPAFB,  Ohio. 

The  effort  by  Sarwan  S.  Sandhu  (P.I.),  his  students  and  Dr.  J.  P.  Fellner  on  the 
research/development  activity  has  resulted  in  the  following  presentations  and  proceedings 
publications: 

“Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling:  Open-Circuit  Voltage 
Equation”  in  the  Proceedings  of  the  8th  International  Conference  on  Electrical  and 
Electronic  Products,  Vol.  34  (2002),  pp.  233-253,  Product  Safety  Corporation, 
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Sissonville,  West  Virginia;  The  Conference  held  at  the  Greenbrier,  White  Sulphur 
Springs,  West  Virginia;  January  14-16,  2002. 

Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling:  Reversible  Open 
Circuit  Voltage  Equation  and  Species  Flux  Equations”,  presented  at  the  8th  International 
Symposium  on  Polymer  Electrolytes”  held  at  the  Eldorado  Hotel,  Santa  Fe,  New  Mexico, 
May  19-24,  2002. 

“Direct  Methanol  Polymer  Electrolyte  Fuel  Cell  Modeling,”  presented  at  the  2002 
AIChE  North  Central  Regional  Student  Conference.  “A  Climate  Change,”  Chem-E-Car 
Competition,  University  of  Michigan,  February  7-9,  2002. 


3  Methods,  Assumptions,  and  Procedures 

Presented  in  this  report  volume,  the  P.I.  has  developed  the  mathematical  equations  to 
predict  the  ionic  current  density  and  the  species  (H+,CH30H,  H2O)  molar  fluxes  for  their 
transport  through  a  solid  polymer  electrolyte  or  polymer-ceramic  material  composite 
membrane  of  a  DMFC  at  the  steady  state,  isothermal  conditions  by  the  application  of 
fundamentals  of  transport  phenomena.  The  developed  transport  equations  are  presented 
in  three  sets.  In  the  transport  set  I  section,  the  equations  are  based  on  the  extension  of 
moderately  dilute  solution  theory.  These  equations  involve  the  Fick’s  law  type  mass 
diffusivities  that  are  known  to  be  strongly  dependent  on  the  species  concentrations.  The 
transport  equation  sets  II  and  III  sections  are  valid  for  very  strong  solutions.  The 
transport  equation  set  II  section  is  based  on  the  application  of  the  generalized  Maxwell- 
Stefan  equations.  The  transport  equation  set  III  section  is  based  on  the  application  of  the 
Lars  Onsager’s  irreversible  thermodynamic  approach  to  transport  processes.  The 
developed  equations  account  for  the  effect  of  the  voltage,  pressure  and  the  species 
concentration  gradients  on  their  transport  fluxes.  These  equations,  describing  species 
transport  fluxes  through  a  solid  polymer  electrolyte  membrane  are  to  be  coupled  with  the 
equations,  yet  to  be  developed,  describing  mass  transfer  through  the  electrode  porous 
backing  layers  and  species  mass  transport  with  electrochemical/chemical  kinetics  in  the 
porous  electrode  layers.  The  coupled  equations,  describing  the  various  phenomena 
occurring  in  a  DMFC  in  operation,  would  be  of  immense  significance  in  the  evaluation  of 
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performance  and  design/development  of  a  DMFC.  The  usefulness  of  the  developed 
transport  equations  may  also  be  seen  in  that  they  can  be  employed  to  experimentally 
determine  the  transport  properties  of  a  membrane  such  as  permeability  of  the  fluid 
mixture  and  effective  mass  diffusivities  of  the  species  by  designing  appropriate  fuel  cell 
experiments. 

4  Results  and  Discussion 

4.1  Flux  Equations  for  a  Direct  Methanol  Fuel  Cell  Solid  Polymer 
Electrolyte  Membrane 

Derivation  of  the  transport  flux  equations  for  a  direct  methanol  fuel  cell  solid  polymer 
electrolyte  membrane  using  i)  moderately  dilute  solution  theory  based  on  Fickian 
diffusion,  ii)  concentrated  solution  theory  based  on  the  generalized  Stefan-Maxwell 
equations,  and  iii)  concentrated  solution  theory  based  on  the  Onsager  irreversible 
thennodynamics  approach  to  transport  processes  are  given  below. 


4.1.1  Transport  Set  I  -  Flux  Equations  using  Moderately  Dilute  Solution 
Theory 

4.1. 1.1  Introduction 

The  transport  equations  given  below  are  quite  general  in  that  they  can  be  applied  to  any 
type  of  membrane  (e.g.  acid  doped  polybenzimidazole,  Nafion  ®  -  DuPont,  a  ceramic 
(clay:  Al2O3.2SiO9.2H2O)-  Nafion®  or  PVDF  (  polyvinylidene  floride)  composite)  with 
some  adjustments  if  required.  The  developed  equations  are  presented  in  the  form  such 
that  transport  and  solubility/permeability  data  either  can  be  obtained  from  experimental 
work  or  available  literature.  It  is  assumed  here  that  the  geometrical  configuration  of  the 
membrane  is  invariant  during  the  isothermal,  steady  state  operation  of  DMFC. 
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4.1. 1.2  Mathematical  Derivation 


In  general,  transport  of  H+  through  membrane  can  be  considered  via  diffusion,  migration 
due  to  electric  field  and  convection  dependent  upon  the  pressure  difference  across  the 
membrane  and  permeability  of  the  membrane  to  the  fluid  mixture  of  CH3OH,  H2O,  H+, 
etc.,  noting  that  the  membrane  penneability  is  function  of  its  molecular  structure  and  void 
fraction,  e  ,  defined  as  the  ratio  of  the  volume  of  the  interstices  between  the  atoms  / 
functional  groups  and  nanometer  /micrometer  pores  to  the  total  membrane  occupied 
volume. 

If  the  gradient  of  the  H+  electrochemical  potential  is  used  as  the  driving  force  for 
diffusion  through  membrane  in  the  presence  of  water  or  water-methanol  mixture,  the  flux 
of  H+  per  unit  area  of  the  membrane  perpendicular  to  the  direction  of  transport  of  H+  can 
be  expressed  as: 

cff+[-^v]  +  c^  (1) 

where  u[!+  =  mobility  of  H+  in  water  or  in  water-methanol  mixture  present  in  the 
membrane,  [(m  /s)/  (J/mol)]; 

t  =  tortuosity  factor  to  account  for  the  zigzag  path  ways  for  the  species  transport 
through  the  membrane; 

e  =  void  fraction  of  interstices,  nanometer/micrometer  size  pores  acting  as 
pathways  for  the  penneation  of  H+,  water,  and  methanol,  if  methanol  crossover  is 
occurring  through  the  membrane; 

cH,  =  fluid  phase  concentration  of  H  ,  [  mol/mf  ]; 

juH+  =  +  ZH+  F®  =  electrochemical  potential  of  H+  in  the  fluid  phase,  [ 

J/mol]; 

v0=  superficial  fluid  phase  velocity,  [m/s]; 

=  chemical  potential  of  H+  in  the  fluid  phase  in  the  presence  of  membrane 
material  [J/mol]; 

zH+  =  number  of  protonic  charges  on  H+=  + 1 ; 
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F  =  Faraday’s  constant,  e.g.,  charge  on  1  g-mol  of  charged  H+  ions  =  96487 
[coulomb/mol]; 

®  =  electrical  potential^  V  ]. 


In  the  first  tenn  on  the  right  hand  side  of  Eq.  (1), 


^  'lU  4  ^ 
H 4 

V  r  J 


is  the  effective  mobility  of  H+. 


This  multiplied  with  (-  V ju  +  j  results  in  the  effective  velocity  for  diffusion  and 


migration  via  electric  field,  and  multiplication  by  the  concentration  c  +  results  in  the 

contribution  to  the  net  flux  NH+ .  The  second  term  on  the  right  hand  side  is  the 
contribution  to  the  flux  via  convective  motion  of  the  fluid  phase. 


The  Nernst-Einstein  relation  is  used  to  express  u n ,  in  tenns  of  D n  [diffusivity  of  H+], 

D  . 

(2) 


H  RT 

where,  D  +  =  H+  ion  diffusivity  in  the  fluid  phase  (in  water  or  water-methanol 

mixture)  in  the  membrane,  [m  /s]; 

T  =  temperature,  [K]; 

R  =  universal  gas  constant,  [8.314 J lmol-K\. 

Replacing  uff+ ,  in  terms  of  D ^  ,  in  Eq.  (1)  leads  to 
eD, 


N, 


H * 


N 


H*  =  D 


H\eff 


(c  \ 

l  RT) 

r^. 

C  A 

^  1 

RT  J[ 

-V^J 

+  C  V 
H *  o 


+  C  .  V 

H *  o 


(3) 

(3a) 


where  D 


D, 


H "  _ 


H+,eff 
2/ 


=  effective  diffusivity  of  H  in  the  fluid  phase  present  in 


the  membrane,  [m  /s]. 


If  the  hydrogen  ion  is  chosen  to  define  the  electric  potential,  O  ; 
Then,  /uH+  =  RT  In 


V  Co  / 


T  z  +FO 

H+ 


(4) 
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The  gradient  of  the  electrochemical  potential  of  any  other  ionic  species  i  can  be  written 
as: 


V//,  =  V 


(5) 


Note  that  in  the  first  term  on  the  right  hand  side  of  Eq.  (5)  the  expression  within  the 
parentheses  is  independent  of  electric  potential  according  to  the  reference  [J.S.  Newman, 
Electrochemical  systems,  Prentice  Hall,  1991]. 


u  -  -L-  Li u  =  RTln  A  - RTIn  A 

z  H  z ...  h 


where  At  is  the  absolute  activity  of  species  i;  zt ,  the  number  of  unitary  charges  on 
species  i. 

a  (  C. A  „ 


A,  =  af  f,  (6a) 

where  af  =  a  proportionality  constant,  independent  of  composition  and  electrical 

state;  function  of  the  nature  of  solute  and  solvent  and  dependent  on  temperature  and 
pressure;  for  condensed  phases  the  pressure  dependence  is  ignored; 
c,  =  molar  concentration  of  species  i,  [mol/liter]; 

c0=  1  mol/liter;  is  introduced  here  to  make  af ,  [a  dimensionless  quantity]; 
f  =  activity  coefficient  of  species,  i,  to  account  for  non-ideal  behavior  of 
species  i  in  the  fluid  phase  present  in  the  membrane,  [dimensionless]. 

(c  \ 

Similarly,  AH+  =  a^+  —  fH+  (6b) 

V  co 

Substituting  for  At  and  A  +  from  Eqs.  (6a)  and  (6b)  into  (6),  one  obtains: 


/6-  —  /V 
V 


& 

H 

In 

f  &  ci 

ai  —Ji 

Zi  In 

l  J 

ZH+ 

«'  c  '  I, 


c,  z. 


RT  (In af - Mn  a„+ )  +  RT  ln^-^ln^  + 


c  z  + 

o  H* 
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RT 


z. 


In/. - —In/ 


(6c) 


v  y 

In  general  for  any  ionic  species,  i,  the  flux  equation  can  be  written  similar  to  Eq.  (3a). 


N, 


f  C;  ^ 


D 


(-Vft) 


+  C  V 

l  o 


kRTj 

V//(  using  Eqs.  (4),  (5),  and  (6c)  is: 


(7) 


V//,  =  V 


RT( In af - Mna/  )  +  RT 


'  c  z  c  '  ' 
In  — - —  In— ^  |  + 


V  C*  V 


°  J 


RT 


A 

J 


+  — -  VjU 
z 


H* 


=  V 


RT  ( In  af 


H * 


Tna^+  )  +  RT 


In— - —  In— ^ 


J 

\  \ 
+ 


V  C°  V  Co  J 


f 


RT 


2; 


A 


In/ - Mn  fH 

V  V  J 


H+ 


(7a) 


+  ^~  V 


f 


-H+  V 


7?7Tn— ^ 

c„ 


+zH+m 


v//,  =  v 


RT(  In  af - Una*  ) 


V 


f  ,c.l 

+  V 

RT  In  — 

l  CJ 

+v 


f 

f  z  Y| 

RT 

In/  -  '  ln/„. 

V 

1  a+  )) 

+  z..F  V  ® 


(7b) 


For  the  condition  of  constant  temperature,  Eq.  (7b)  becomes 


V/u,  =  RT  V  In  ( —  )  +  ztF  V  ®  +  RT  V 
c„ 


In/ - —In/ 


(7c) 


Then,  the  flux  Eq.  (7)  for  the  constant  temperature,  can  be  written  as 

7,  =  -A.«r  U,V<I>-  D,,g  c, V ln(  — ) 


"  A',e/f  C/^7 


In/ - —  hi/^ 

V  y 


+  C-  V 
*  o 


(8) 
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Simplification  of  the  second  term  on  the  right  hand  side  of  the  Eq.  (8)  leads  to 


r  zF 

N .  =  -  D  — — c  VO  —  D  ff  Vc 

»  ’,eff  1  ’> eff  ' 


"  A',<#  Ci ^ 


Z.. 


A 


In  / - -In  fH 

ZH+  J 


+  C,  V„ 


(8a) 


In/.  -- i-In/^  =  In/  -  In/; 


(8b) 


In/ - —  ln/„+  =  In 


H+ 


(  \ 
f, 


f  »- 

\J  H*  J 


Using  the  above  result  in  Eq.  (8a),  the  following  equation  is  obtained. 
zF 


Ni  =-  Di,eff  C,VO-  Di,eff  VU 


Di,eff  c,Vln  f.H+  +  civ0 


where  / 


r  \ 
ft 


f  H+ 

\jh +  y 


(8c) 


(8d) 


Either  from  Eqs.  (3a)  and  (4)  or  from  Eq.  (8d)  directly  with  i  =  H+,  one  obtains  the  flux 


equation  for  H+  ions  : 

T  Zh+F  -  - 

N  +  =-  D  +  - c  +VO—  D  Vc  +  +  c  +  v 

H+  H  ,eff  H*  H*  ,eff  H *  H*  » 


(9) 


Note  that  the  activity  coefficient  /  of  the  charged  species  can  be  estimated  following  the 
procedure  given  in  the  literature  (Ref.  J.S.  Newman,  Electrochemical  systems,  pp  86-115, 
Prentice  Hall,  1991  edition). 


The  superficial  velocity  vo ,  in  Eqs.  (8d)  and  (9),  is  the  volume  rate  of  flow  through  a 

unit  cross-sectional  area  of  the  solid  membrane  material  plus  fluid  phase,  here,  composed 
of  water,  methanol  and  H  .  It  is  averaged  over  a  small  region  of  the  membrane  facial  area 
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perpendicular  to  the  flow  but  large  with  regard  to  the  interstitial  micro  size  pores.  The 
velocity  va  can  be  described  by: 

Vp-pg-p  — -  (10) 

p 

V  )_ 

In  this  equation,  //  is  the  fluid  viscosity  [kg/ms];  kpjM ,  the  fluid  phase  permeability 

through  the  membrane,  \m2  J .  The  first  term  on  the  right  side  of  Eq.  (10)  is  the  pressure 
gradient  related  force  on  the  fluid  phase,  the  gravitational  force  is  the  second  term,  and 
the  third  term  is  the  net  force  per  unit  fluid  phase  volume  due  to  electric  field  effect.  This 
tenn  may  be  labeled  as  the  driving  force  for  the  electro-osmotic  transport  of  the  fluid 
phase  through  the  membrane.  Also  note  that  c,  is  the  molar  concentration  of  the  ionic 

species,  i  (with  the  charge  number  z. )  in  the  fluid  phase,  with  p  being  the  fluid  phase 
mass  density. 


For  the  transport  of  the  fluid  phase  perpendicular  to  the  direction  of  the  gravitational 
field,  or,  in  the  situation  of  the  second  tenn  on  the  right  hand  side  of  Eq.  (10)  being 
negligibly  small  relative  to  the  first  and  the  third  terms;  Eq.  (10)  is  reduced  to: 

vo=  .  hlL  yp-2>,F(-V®)c,  (11) 

M  \_  i 

In  the  event,  either  (- Vo)  electric  field  is  negligibly  small  in  the  membrane  space,  or 
there  is  no  electric  field  in  action,  then  Eq.  (11)  becomes 

v„=-— M  (1  la) 

Note  that  this  is  the  usual  Darcy’s  law  for  the  flow  through  a  porous  medium  with  the 
pressure  gradient  as  the  only  driving  force  for  the  fluid  movement.  In  the  case  of  the  fuel 
cell,  if  Vp  is  large  because  of  the  difference  in  the  total  pressures  in  the  anode  and  the 
cathode  compartments  and  the  small  thickness  of  the  membrane,  va  may  acquire  a 
relatively  significant  value  depending  upon  the  value  of  kpfM  and  //  of  the  fluid  phase 
mixture  of  water,  methanol,  and  H+  .  In  the  event,  Vp  is  negligible  either  due  to  equal 
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total  pressures  in  the  fuel  cell  anodic  and  cathodic  compartments  or  due  to  a  relatively 
thick  membrane,  then  Eq.  (11)  becomes: 


(11a) 


(Hb) 


if  it  is  assumed  that  the  fluid  phase  is  electrically  neutral,  then, 

I  MV  =» 


(lie) 


Consequently,  Eq.  (lib)  leads  to:  vo  =  0  ,  i.e.,  there  is  no  effect  of  the  electric  field  on  the 

fluid  convective  motion  through  the  membrane.  This  is  not  the  valid  conclusion  because 
the  literature  does  give  evidence  of  the  fluid  transport  through  the  membrane  under  the 
field  effect.  Then,  we  have  to  accept  that  the  membrane  material  to  be  charge  neutral,  but 
the  fluid  phase  permeating  through  membrane  phase  is  not  charge  neutral.  The  net 
charge  per  unit  fluid  phase  volume  is  ( zH+cH+F  ).  Now,  Eq.  (lib)  is  written  as: 

V„=  [k— )  (-vo)  [zh.ch.f)  (lid) 

Note  that  vo  may  have  relatively  significant  value  depending  on  the  values  of  (-  Vo) 
[electric  field  strength  in  the  fluid  phase  in  the  membrane],  c  +  in  the  fluid  phase  in  the 
membrane  and  k pfM  . 


For  the  charge  transport  by  the  only  H+  ions  in  the  DMFC;  Eq.  (11)  can  be  written  as 

V„=  -f— <12> 

V  F  ) 


Substituting  for  vo  from  Eq.  (12)  into  Eq.  (9)  and  simplification  leads  to 


—W  (-  Vp)  (13) 

\  M  ) 
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Note  the  structure  of  the  individual  tenns  contributing  to  the  flux  of  H+  ion  through  the 
membrane.  In  the  first  term  on  the  right  hand  side  of  the  Eq.  (13),  (-  Vo),  the  electric 
field  is  acting  as  the  driving  force  and  (-  Vc  +  )  and  (-  Vpj  are  the  concentration  and 

pressure  gradient  driving  forces  in  the  second  and  third  terms,  respectively. 

Substitute  for  vo  from  Eq.  (12)  into  Eq.  (8d)  to  obtain: 


This  equation  is  valid  for  any  ionic  and  non-ionic  species. 


Methanol  and  water  are  charge  neutral  species,  i.e.,  zCHOH  =  zH  0  =  0 ,  i.e.,  number  of 
charges  on  each  of  these  species  is  zero. 


Also,  / 


/  (activity  coefficients  of  the  charged  neutral  species  i) 


Equation  (14)  is  written  for  a  charge  neutral  species  i  as: 

-  [kjf\z»-Fc«-c'  (-^®)+  (-^)  + 

\  M  J 


Di,eff  Ci  I” 


,  (- VlnX) 


For  the  situation  where  it  is  reasonable  to  assume  that  (v  In/)  =  0,  i.e.,  variation  of  f 
across  the  membrane  is  negligible;  Eq.  (15)  becomes: 

DIJff  (-Vc,)+f^)c,(-Vp)  (16) 

V  A  7  v  ^  J 
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where,  i  =  CH3OH,  HoO. 


Note  that  the  first  term  on  the  right  hand  side  of  Eq.  (16)  involves  (-  VO  J,  electric  field 
as  the  driving  force  even  for  a  charge  neutral  species  such  as  CH3OH,  H2O,  for  example. 
Here  this  term  is  labeled  as  the  electro-osmotic  contribution  to  the  flux  of  charge  neutral 
species  i  through  the  membrane.  The  second  tenn  on  the  right  side  is  the  diffusion 
contribution  to  the  flux  of  species  i  and  the  third  term  accounts  for  the  flux  contribution 
by  the  fluid  motion  set  up  by  (-  Vp),  the  pressure  gradient  in  the  case  of  different  total 
pressures  in  the  anode  and  cathode  compartments.  Note  that  the  electro-osmotic  flux 
contribution  does  depend  on  the  fluid  phase  viscosity,  p  ,  its  permeability  through 
membrane,  kpi\i,  and  H+  ion  molar  concentration  in  the  fluid  phase  in  the  membrane 
matrix. 


4.1. 1.3  One-Dimensional,  Steady-State,  and  Isothermal 

We  shall  now  concentrate  on  obtaining  the  equations  for  the  transport  of  species  H+, 
CH3OH,  and  H2O  through  the  membrane  under  the  steady  state,  isothennal  conditions  for 
the  one-dimensional  (x-)  case. 


The  x-components  of  Eqs.  (13)  and  (16)  are: 


Zh*F 
v  RT  J 


FH*,effCH *  + 


( k  \ 

pJM 

l  F  2 


z  Fc 


dQ> 

V  dxj 


+  D, 


H\eff 


dc. 


v  j 


+ 


f  kpJM  "] 

c 

(  dp ^ 

l  F  j 

lh+ 

l 

N, ..  = 


f  kpjM  (  30^ 

zh*Fch*c 


v  F  J 


V 


+  D 


i,eff 


f  dc .  ^ 
V 


+ 


f  kpJM  ' 

F  J 


/  dp ^ 
V  dxJ 


(17a) 

(17b) 


(molar  flux  for  a  charge  neutral  species,  i  =  CH3OH,  H20) 

The  current  density  per  unit  area,  perpendicular  to  x-coordinate  in  the  membrane,  is 
given  by: 
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1=1 


j,  =  £  bNU,]  (18) 

(18a) 
(18b) 

Note  that  the  H+  ion  is  the  only  charge  carrier  through  the  membrane  so  that  the 
transference  number  of  the  species  H+  is  one. 

Using  Eq.  (17a)  for  N H+  x  in  Eq.  (18b), 


=  z  F  N  +  z  FN  +  z  FN 

H*  H+,x  ^CH^OH1  JV  CH3OH,x  ^ H-,01  ly  H20,x 


=  ZH^NH\x 


Jx 


( z1  F2^ 

AH+r 

RT 


DH\effCH+  + 


KpfM 

v  M  j 


z2  +  F2c2  + 

H*  H  + 


'  d®' 

V  dx  J 


+  ZH^DH+,eff 


dc. 


dx 


+ 


fir  ^ 

KpjM  1 

f*  7  F 

(  dp ^ 

l  M  j 

ch+  zm* 

l  dx) 

(19) 


( z2  F2^ 

ZH+r 

RT 


D 


H* 


V  dxj 


+ 


1 pm 

v  M  j 


2/72  2 
VF  CH 


+ 


ZH^DH\eff 


dc. 


dx 


+ 


fk  ' 

lpJM 


CH+  VF 


V  dx  j 

/  dp2 
V  dx) 


(20) 


v  —  j  \  F  j 

Equation  (20)  requires  the  explanation  to  various  tenns.  The  first  term  on  the  right  hand 
side  expresses  transport  via  H+  ion  movement  due  to  the  action  of  the  electric  field  in  the 


x-direction.  The  coefficient  of 


f  fclU 


V  dx  j 

through  fluid  phase  present  in  the  membrane  matrix,  i.e., 


is  here  defined  as  the  electric  conductance  of  H 


2  /-»  2 

Kelec  =  VF 


D. 


H\eff 

RT 


CH+  =  ZVFl  V 


(20a) 


where  u 


H + 


ion  mobility  in  the  fluid  composed  of 

{ciI  .OIT  H20,  and  H  ^present  in  the  membrane 


V  ^  J 


(20b) 


The  second  term  on  the  right-hand  side  of  the  Eq.  (20)  explains  the  current  density 
contribution  associated  with  the  fluid  “bulk  motion”  created  by  the  action  of  electric  field 
on  H+  ions  present  in  the  fluid.  Third  term  is  the  current  density  contribution  because  of 
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gradient  in  the  H+  ion  concentration.  Fourth  term  is  the  contribution  to  the  current  density 
if  the  pressure  gradient  prevails  for  the  different  total  pressures  across  the  membrane  in 
the  electrode  chambers. 

Equation  (20)  may  now  be  written  as 


J  x  = 


K 


elec 


'  dO  ' 

+ 

j  k  pJM  'j 

(z  Fc  V 

(  dO'' 

l  J 

l  M  J 

|  V  H+  H+  > 

l  dx  J 

z  +FD  + 

H +  H  ,eff 


dc 


H  + 


dx 


+ 


f  u  \ 
lP.fM 

v  M  , 


+ 


dp 

dx 


(21) 


For  the  medium  of  assumed  uniform  dielectric  constant  (here  that  of  water  or  of  water- 
methanol  mixture  present  in  the  membrane),  the  Poisson  equation  is  used  to  relate  the 
Laplacian  of  the  electric  potential  to  the  charge  density  in  the  fluid  phase  in  the 
membrane: 


F 


V2®  = - Yj 


F 


twl’ 


where  e  =  pennittivity  of  the  medium  (here  that  of  the  fluid  phase) 
=  (dielectric  constant)  (e0 ,  permittivity  of  the  free  space) 

f  coulomb  ' 


=  (dielectric  constant) 


8.8542  *  10 


-14 


V 


=  (dielectric  constant  of  the  fluid  phase) 
F  =  Faraday’s  constant  =  96487 


f  coulomb ^ 


volt  -  cm  j 
8.8542  *  10~12 
coulomb 


V 


volt  -  m 


mole  of  unitary  charges 
For  the  one-dimensional  (x-)  case,  Eq.  (22a)  becomes: 

J2®  F 


dx1 


[w. 


(22) 

(22a) 


(22b) 


The  material  balance  or  continuity  equation  for  a  component  i  for  the  steady  state 
condition,  with  no  reaction,  is: 
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V  •  TV,.  =  0 


(23) 


dNix  dNiy  dNt 


=  0 

dx  dy  dz 

For  the  one-dimensional  component  in  the  x-direction, 

dN , 


(23a) 


l,X  _ 


dx 


=  0 


(23b) 


Using  Eq.  (17a),  the  following  equation  is  obtained  from  Eq.  (23b)  under  the  assumption 

f  dO  h  f  ^ 


of  constant 


V  dx  j 


dp 

\dx  j 


,  kpfM ,  //  and  D/r  for  the  isothennal  cell  operation. 


Zh+F 

V  ^  , 


D. 


dc , 


H  ~eff  dx 


■  + 


fir  \ 
pjM 

P  J 


dCH  + 

Z»-F2C»-^k 


dO 

dx 


+  D. 


H\eff 


d2c 


H* 


dx2 


+ 


r  kPfM" 

dCH+ 

(  _dp^ 

l  M  ) 

dx 

v  dx  j 

0 


(23c) 


The  modified  continuity  equation  for  the  fluid  mixture  of  (H20,  CH3OH,  and  H+)  through 
the  nano/micro  porous  medium  is  applied  for  the  permeation  of  fluid  phase  mixture 
through  the  electrolyte  membrane  at  the  steady  state  conditions. 

V*bvo]=0,  (24) 

where  p  =  fluid-phase  density. 

For  the  one-dimensional  (x-component)  case,  Eq.  (24)  is  written  as: 
d 


dx 


(u„J  =  0 


(24a) 


The  solution  of  Eq.  (24a)  is 
(pvo.x  L=0  =(pvo,X  ■  =  m„  =  constant, 


i,x=Su 


(24b) 


where  ma  =  superficial  fluid  mass  flux  permeating  through  the  membrane  along 

x-axis  per  unit  time  per  unit  area  of  the  membrane  perpendicular  to  the 
x-axis. 

The  x-component  of  Eq.  (12)  is 


v0>*  = 


fir  A 

pjM 


J  L 


dp 

dx 


+  VFV 


dO 

dx 


(24c) 
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The  fluid  phase  mass  flux  is  related  to  the  pressure  and  electrical  gradients  as  follows: 


=  P  v0 . 


pk 


pm 

P  j 


dp 

V  dx  j 


+  VFV 


dd> 

dx 


(24d) 


For  the  electric  potential  distribution,  Eq.  (22b)  is  written  as 


d2P 

dx2 


F 


twl’ 


(25) 


where  c  +  is  the  H+  ion  average  concentration  defined  as 


x=5u 


JcH+  dx 


c  =  x=0 


S 


(25a) 


M 


where  SM  =  thickness  of  the  membrane.  In  fact,  one  may  approximate  this  quite 
reasonably  by 


c 


h * 


+  c , 


(26) 


where  c  2r  is  the  H+  ion  concentration  in  the  fluid  phase  at  the  anode-membrane  interface 
and  c'e  the  concentration  of  H+  at  the  membrane-cathode  interface  on  the  membrane 

side.  For  the  situation  of  relatively  fast  electrode  kinetics  associated  with  the 
electrochemical  reactions,  especially  fast  kinetics  at  the  cathode,  the  H+  ion  concentration 
at  the  membrane-cathode  interface,  c(ir  may  be  approximately  set  equal  to  zero.  Thus, 


c 


H+ 


c 


A 

H+ 


2 


Equation  (25)  is  solved  using  the  following  boundary  conditions: 
At  x  =  0,  ®  =  cf> 4  (electric  potential  in  the  fluid  phase  at  the 
anode-membrane  interface) 

At  x  =  SM ,  ®  =  ®  c  (electric  potential  in  the  fluid  phase  at  the 

cathode-membrane  interface) 

The  solution  to  Eq.  (25)  is 


f  \ 


(O-®^1)  =(®c  -®^) 


\^M  J 


+ 


z  +c  F 

H*  H* 

2  e 


(SMx-x2) 


(26a) 


(27a) 

(27b) 


(28) 
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Note:  Here,  ( 0  -  <D  )  =  -(absolute  value  of  the  cell  voltage) 

=  -|oc  -  O'4 1 

Note  that  Eq.  (28)  describes  the  electric  potential  in  the  fluid  present  in  the  interstices  in 
the  membrane  matrix. 

From  Eq.  (28)  one  obtains 

0  c  -  0 


J0 

dx 


V  J 


+ 


Z~^~  fe-2x) 

v  2e 


(29) 


0C  -0 


5 


+ 


M 


v  2e  y 


{SM  ~  2x) 


(29a) 


The  result  in  Eq.  (29)  is  used  in  Eq.  (24d)  and  the  resultant  equation  is  solved  for  the 
pressure  distribution  in  the  fluid  phase  in  the  membrane  matrix  to  obtain 

( P-P'A)= 


m0n8i 


M 


V  kpjMP  j 


f  \ 


v  8 m  j 


—  7  c  F 

ZH*'H+r 


(0C  -0"*) 


f  -  \  ( z  c  FS2  ^ 

ZH+LH+rUM 


\8m  j 


+ 


2  e 


f  ..  \ 


2  A 


8 


M 


\8 M  J 


(30) 


where  ( -  0^ )  =  -(absolute  value  of  the  cell  voltage), 


t.A 


=  total  pressure  at  the  anode -membrane  interface, 


A  C 

p  =  average  fluid  density  = - ,  where  p  A  and  p  are  the  fluid  densities  at 


the  anode -membrane  and  cathode-membrane  interfaces,  respectively. 

Using  the  boundary  condition  of  p  =  p  ’  (total  pressure  at  the  membrane-cathode 
interface)  at  x  =  SM,  one  obtains  the  following  result  from  Eq.  (30): 


m„  = 


F  ) 


pLA-p,c 

8„ 


0C  -  04 


+  WF- 


8 


M 


(31) 


This  equation  shows  the  effect  of  the  pressure  and  electric  potential  differences  on  the 

penneation  flux  of  the  fluid  phase  through  the  membrane.  For  the  case  of  pL  A  (total 

t  c 

pressure  at  the  anode)  greater  than  p  ’  (total  pressure  at  the  cathode),  both  of  these 
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driving  forces  are  seen  to  cause  the  fluid  phase  to  move  towards  the  cathode.  For  the 

case  p1'  A<  p1' c,  the  effect  of  these  driving  forces  on  m0  is  opposite.  For  the  case  of  p1,  A  = 

t  c 

p’  ,  Eq.  (31)  becomes: 

(32) 

This  simple  result  tells  us  that  due  to  the  electric  field  effect  there  is  a  flux  of  fluid 
mixture  that  permeates  through  the  membrane.  For  a  membrane  of  given  molecular 
configuration,  one  way  to  decrease  ma ;  therefore,  to  decrease  the  flux  of  methanol  is  to 

increase  the  membrane  thickness.  Other  way  is  to  use  a  membrane  of  such  a  molecular 
structure/configuration  that  the  penneability,  kpfM,  is  decreased  in  such  a  way  that 
penneation  of  methanol  is  selectively  reduced  more  than  water  with  either  no  or 
enhancement  effect  on  the  penneation  of  H+  ions. 

Using  Eq.  (15),  one-dimensional  molar  fluxes  of  H20  and  CH3OH  through  the  membrane 
are 


N, 


CH,OH,x 


V  ) 


'  CH2OH  '  H+ 


dO 

dx 


+  D. 


CH,OH,eff 


dc 


CHxOH 


dx 


+ 


(electro-osmotic  transport) 


(diffusion  transport) 


+ 


( k  \ 

KpJM 

(  _dp^ 

{  V  ) 

cch3oh  | 

(  dx  j 

(33) 


(pressure  gradient  transport) 
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f  dO  ' 

+  n 

c 

dCH20 

l  V  ) 

CH2Q\Z H*tCH*  ' 

v  dx  j 

^ H20,eff 

dx 
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(electro-osmotic  transport) 


(diffusion  transport) 
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V  ) 


'H20 


f  dp A 
V  dx  j 


(34) 
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(pressure  gradient  transport) 


For  the  condition  of  the  steady  state,  the  local  continuity  equations  for  H2O  and  CH3OH 
are 


( dN  ^ 

CH,OH,x 

dx  j 


=  0 


(35a) 


v  dx  j 


0 


(35b) 


Substitution  for  NCHOHx  and  N n  o  x  from  Eqs.  (33)  and  (34)  into  Eqs.  (35a)  and  (35b), 
respectively,  and  simplification  leads  to 


(k  > 

^  pfM 

d  " 

l  F  J 

dx 

CCHiOH  (Z h*  F C H +  ) 
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K  pJM 

d 

( 
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l  »  J 

dx 

uCH3OH  1 

dx  J 

=  0 


(36a) 
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K pfM 

d  " 

l  »  J 

dx 
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f  cMF 

y  dx  j 


+  D 


6  d2c  A 

u  lh,o 


H20,eff 


dx  j 
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fk  \ 

K  pfM 

d 

r> 

r  dP  \ 

l  F  ) 

dx 

uH20 

y  dx  J 

(36b) 


From  the  Eq.  (24d),  with  the  assumption  of  constant 


l  F  J 


and  constant  fluid  phase 


density  p  in  the  membrane  phase,  the  following  equations  are  obtained: 
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f  k  3 

K  pjM 
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f  ■  \ 

ma 

dcH2o 
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d 

/> 

dp\ 
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(37b) 


Substitutions  are  made  from  Eqs.  (37a)  and  (37b)  for 


r  k  3 

K  pjM 
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dx 

'CH,OH\^h 
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d<t> 

dx 


f . 


and 


A 
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v  P  J 


dx 
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H* 


dO 

dx 


into 


Eqs. (36a)  and  (36b),  respectively,  to  obtain,  after  simplification,  the  following  results: 


d2c 


CH,OH 


dx 2 


pD 


dc 


CH,OH  _ 


CH2OH  ,ejf  ) 


dx 


=  0 


(38a) 


dx1 


ni. 


pD 


dc 


H-,0 


H20,eff  J 


dx 


=  0 


(38b) 


These  ordinary  differential  equations  are  solved  using  the  following  boundary  conditions: 
At  x  =  0  (membrane-anode  interface); 

CCH2OH  —  CCH3OH,x=0  ’  CH20  —  CH2Q,x= 0  (39a, b) 
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At  x  =  8m  (membrane-cathode  interface); 


_  M 

CCH3OH  CCH3OH,x=Sm  ’ 


_  M 

CH20  CH20,x=Su 


(40a, b) 


The  solution  of  the  differential  Eqs.  (38a)  and  (38b)  based  on  the  boundary  conditions 
Eqs.  (39)  and  (40),  is 


_  M 

CCH3OH  CCH3OH,x=S3 
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where  ax  = 
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pD , 


CH3OH  ,eff 


(41) 

(41a) 
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CH20  C H20,x=Sm  ^  V  CH20,x= 0  "  CH20,x=8m 
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a2  &M  \ 


1  0a2^M 

V  i  -  e  J 


where  a2  = 


m , 


pD 


H20,eff 


(42) 

(42a) 


Methanol  and  water  concentrations  in  the  membrane  matrix  at  the  interfaces  with  the 
porous  electrodes  are  now  related  to  their  concentrations  in  the  fluid  mixtures  in  the 
electrodes  using  the  phase  equilibrium  condition. 

At  the  membrane-anode  interface,  x  =  0; 


a 


M 

CH3OH,x=0 


a 


A 

CH3OH,x=0  > 


(43a) 


a 


M 

H20,x= 0 


=  a 


A 

H20,x= 0  ’ 


(43b) 


where  a;  =  activity  of  species  i. 

Expressing  the  activities  in  terms  of  molar  concentrations, 
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(44a) 


M  rM 

cch3oh,x=oJ  CH3OH 

~cA  fA 

uCH3OH,x=OJ  CH3OH 

Co 

c° 

M  rM 

ch2o,x=oJ  h2o 

A  rA 

ch2o,x=oJ h2o 

Co 

Co 

(44b) 


where  c,  =  molar  concentration  of  a  species, 

fM  =  activity  coefficient  of  species  i  in  the  fluid  phase  in  the 
membrane  phase, 

fA  =  activity  coefficient  of  species  i  in  the  fluid  phase  in  the  porous 
anode, 

c0  =  nonnalization  concentration  =  1  (mol/liter). 


Equations  (44a)  and  (44b)  leads  to 
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M 

CH3OH,x= 0 
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J  CH3OH 
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CH,OH  j 
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M 

H20,x= 0 


(  r  A  \ 
A  JH20 

CH2O,x=0  rM 
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(45b) 


Similarly,  at  the  membrane-cathode  interface  at  x  =  8 M  , 
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J  h2o 

rM 

yJ h2o  j 


At  first  approximation;  particularly,  for  the  situation  of  non-availability  of  the 
in  the  membrane  phase,  one  may  take: 


(46b) 
f]M  data 
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f/.o 

— — 
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— - — 

fM 

\Ji  J 

\J i  y 

=  1. 


From  Eq.  (24d), 
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^  dx  J 

On  substitution  for 
after  simplification. 
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from  Eq.  (47)  into  Eqs.  (33)  and  (34)  to  obtain, 


N, 


M 

CH3OH  ,x " 


^  m.  ^ 


P  J 


c  +  D 

'-'CH2OH  ^  CH2OH  ,ejf 


f  dc  ^ 

CHMH 


dx 


N 


M 

H20,x  ' 


^ m  ^ 


v  P  ) 

From  Eqs.  (41)  and  (42), 
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and 
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dx 


from  Eqs.  (49a)  and  (49b)  into  Eqs.  (48a) 


J 


v  '"~v  j  v 

and  (48b),  respectively,  the  following  expressions  are  obtained  for  methanol  and  water 
fluxes: 
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Note  that  the  fluid  phase  flux  through  the  membrane  matrix  is  given  by  Eq.  (3 1). 

Assuming  the  fluid  phase  permeating  through  the  membrane  matrix  to  be  consisted  of  H+, 
water  and  methanol  (neglecting  dissolved  species  such  as  nitrogen,  oxygen), 


=  1 

i 

where  Xj  is  the  mole  fraction  of  species  i.  Therefore, 

=  1 


(52) 


C  H*  CCH,OH  CH20 


(52a) 


where  c H,  ,  cCH  OH ,  c  HO ,  and  cT  are  the  molar  concentrations  of  H  ,  methanol,  water 

and  fluid  mixture,  respectively.  Here  it  is  assumed  that  the  fluid  phase  total  molar 
concentration  is  constant  throughout  the  membrane  matrix.  Differentiating  Eq.  (52a) 
leads  to 
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,  (53) 

dx  dx 

Using  the  expressions  from  Eqs.  (49a)  and  (49b)  in  Eq.  (53),  the  following  expression  for 

rdc  \ 

— —  is  obtained. 
dx 
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(54) 


v1  "  J 

From  Eq.  (52a), 

(55) 

Applying  this  equation  at  x  =  0  (membrane  -anode  interface)  and  at  x  =  SM  (membrane- 


M  /  M 

CH+  =CT-  \C 


cathode  interface),  the  following  expressions  for  c1^  i  ()  and  c  ]'.  i,  are  obtained. 


" H\x=Su 
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M  _  l  M  .  M  \ 

CH*,x= 0  CT  "  \CCH3OH,x= 0  +CH2O,x=0) 

M  =  r  (rM  ,  M 

°H\x=Sm  CT  ‘  \CCH3OH,x=Su  +CH20,x=Sm 


(56a) 

(56b) 


Then,  the  average  concentration  of  hydrogen  ion  cH+  in  the  membrane  matrix 
(see  Eq.  26)  is  given  as 
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CH+  +CH+ 
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f  M  ,  M  M  M 
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(57b) 


The  methanol  and  water  concentrations  in  the  membrane  at  the  membrane-porous 
electrode  interfaces  are  to  be  computed  from  the  concentrations  of  these  species  in  the 
porous  electrodes  as  shown  in  Eqs.  (45)  and  (46). 

The  hydrogen  ion  concentration,  in  the  membrane  as  a  function  of  distance  can  be 


calculated  using  the  Eqs.  (45),  (46),  (54)  and  (55). 
From  Eqs.  (29a)  and  (30), 
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mjA 

Wj  $m  2e 


Because  of  the  steady-state  condition,  the  following  equations  are  obtained  using  Eqs. 
(17a)  and  (21): 
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(65) 


/77  772 

Note  that  cq  = - 2 - ;  a2  = - 2 - . 

P^CH3OH,eff  PP*H20,eff 

Also,  cIJ+  =  c, M+  is  to  be  obtained  from  Eq.  (57b),  and  (c\Tl  is  to  be  obtained  from 

ti  H  \  H  /@x=0 

Eq.  (56a).  Of  course,  m0  is  to  be  obtained  from  Eq.  (31). 


Equations  (63)  and  (65)  can  be  used  to  compute  the  hydrogen  ion  flux  and  current 
density,  respectively,  through  the  fluid  phase  in  the  membrane  matrix. 


4.1.2  Transport  Set  II  Flux  Equations  using  Maxwell-Stefan  Relations 

4.1. 2.1  Introduction 

An  extremely  wide  variety  of  membrane  types  exist.  One  may  distinguish  between 
biological  and  synthetic  membranes.  Each  type  has  wide  range  of  composition  and 
behavior.  The  synthetic  membranes  may  be  classified  into  (a)  “homogenous”,  (b) 
“microporous”  and  (c)  mixtures  of  these  two  types.  In  a  homogenous  membrane  the 
matrix  plays  the  role  of  a  true  solvent  for  permeating  species.  The  permeating  species  are 
confined  to  matrix-free  regions  in  microporous  membranes.  These  aspects  are  important 
for  materials  viewpoint.  However,  the  formalisms  needed  to  explain  the  transport  of 
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species  through  them  is  much  the  same.  Also  note  that  the  membranes  are  used  in  a  wide 
range  of  process  conditions. 


Anode-membrane 

interface 


Cathode-  membrane 
interface 


Legend: 


x  =  spatial  coordinate 
SM  =  thickness  of  the  membrane 


Figure  1.  Polymer  electrolyte  membrane  separating  the  porous  anode  and 
cathode  layers  (sketch  not  to  scale). 

4.1. 2.2  Mathematical  Derivation 

The  species  transport  equations  given  below  are  quite  general  in  that  they  can  be  applied 
to  any  type  of  membrane  (for  example,  acid  doped  polybenzimidazole,  Nafion  ®  ,  a 
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ceramic  (clay:  A1 2  O  3 .2SiO  2 ,2H  2  O)-  Najlon ®  or  PVDF  (  polyvinylidene  floride) 

composite).  It  is  here  assumed  that  the  geometrical  configuration  of  the  membrane  is 
invariant  during  the  isothermal,  steady-state  operation  of  the  direct  methanol  fuel  cell 
(DMFC). 

Here  a  membrane  is  viewed  as  an  insoluble  selectively  penneable  matrix  M  and  one  or 

more  mobile  permeating  species;  a,/3, . The  membrane  is  defined  here  as  having 

negligible  curvature,  i.e.,  its  thickness  is  much  less  than  its  surface  radius  of  curvature. 
This,  in  turn,  implies  that  the  mass  transport  is  uni-directional  and  perpendicular  to  the 
membrane  surface.  The  membrane  matrix  is  immobile.  That  is,  the  velocity  of  the 
matrix  is  zero  with  regard  to  a  stationary  laboratory  coordinate  system.  Velocities  of  the 
mobile  species  are  referred  to  the  coordinate  system  fixed  in  the  stationary  membrane. 
The  diffusional  response  times  within  the  membrane  are  assumed  to  be  short  in 
comparison  with  those  in  the  mixtures  in  contact  with  its  boundaries  (located  at  x  =  0  and 
x  =  8  as  shown  in  the  sketch;  where  8  is  the  thickness  of  the  membrane).  This 
effectively  means  that  for  every  species  a  pseudosteady  state  prevails  in  the  membrane 
even  when  slow  transient  conditions  prevail  in  the  mixtures  in  contact  with  its 
boundaries. 

The  Maxwell-Stefan  equations  given  below  satisfy  the  above  mentioned  membrane 
constraints  and  are  specialized  to  generate  reliable  and  compact  descriptions  of 
multicomponent  transport  in  the  membrane  located  between  the  electrodes  of  a  direct 
methanol  fuel  cell. 

It  is  recognized  that  the  matrix  is  one  of  the  components  of  the  mixture  of  matrix,  M  and 
other  electrochemical/chemical  species  such  as  H  + ,  CH3OH,  H2  O,  etc.  The  Maxwell- 

Stefan  equations  are  applied  only  for  the  mobile  species  (Ref.  Bird  et  al. [2002]).  From  a 
thermodynamic  viewpoint,  we  are  concerned  with  the  number  of  independent  mobile 
components  in  the  mixture  bathing  the  membrane  because  it  is  the  external  solution  that 
detennines  the  equilibrium  state  of  the  membrane.  As  mentioned  in  the  above  reference, 
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determination  of  the  effective  molecular  weight  of  the  membrane  matrix  is  very  difficult 
for  most  situations,  the  internal  system  is  defined  to  include  only  the  mobile  species  and 
their  mole  fractions  are  defined  to  sum  to  unity.  However,  this  restriction  may  be  relaxed 
if  the  average  molecular  weight  of  the  membrane  material  and  its  actual  membrane 
material,  M,  mass  density  for  its  configuration  between  the  fuel  cell  electrodes  is  known. 
The  interaction  of  each  mobile  species  with  the  membrane  material,  M,  is  quite 
significant;  particularly,  in  the  case  of  a  membrane  with  the  interstices  of  nanometer  scale 
or  less.  Therefore,  the  modified  mass  diffusivity  parameter  associated  with  diffusion  of  a 
species  a  accounting  for  its  interaction  with  the  membrane  is  defined  as 


1 


B, 


B, 


(66) 


«M  «M 

where  xM,  BaM  are  the  membrane  material  mole  fraction  and  mass  diffusivity  of  the 

component  a  through  membrane,  M;  as  it  appears  in  the  Maxwell-Stefan  equations 
without  excluding  the  membrane  material;  respectively.  The  generalized  set  of  the 
Maxwell-Stefan  equations  for  the  fuel  cell  application  is: 


allspecies  v  v 

V  XocX p  -  \ 

L,  -^—\v«-vn) 

P=\  Vap 


■xaVlnaa-  xa 


fva) 

, 

(  F  \ 

Vp-  Xn  z 

l RT ) 

x  (X  u 

krt) 

VO 


(67) 


(Here  a  =  mobile  species  only  =  H  + ,  CH  3  OH,  H  2  O  only) 

Note  that  according  to  the  procedures  used  in  the  reference  mentioned  above,  equation 
(67)  should  be  used  for  the  species  that  are  mobile  and  penneate  through  the  membrane. 
However,  sum  on  the  left-side  of  Eq.  (67)  should  be  carried  over  all  species  including  the 
membrane  matrix,  M.  In  the  event  of  that  mole  fraction  of  the  membrane  matrix  M 
cannot  be  detennined  with  a  reasonable  accuracy,  it  has  been  suggested,  in  the  reference 

X  1 

given  above  that  -  —  be  replaced  by  — T —  in  Eq.(67)  and  the  species  mole  fractions  be 
based  on  the  mixture  molar  concentrations  excluding  the  membrane  matrix  material. 


In  this  equation,  xa ,  x p  are  the  mole  fractions  of  the  species  a  and  ft  ;  va ,  vf, ,  the 
resultant  velocities  of  the  species  a  and  /?  ,  respectively  ;  Bap ,  the  diffusion  coefficient 
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for  the  interaction  of  species  a  with  the  species  jB  ;  va ,  the  partial  molar  volume  of 
species  a  in  the  membrane  ;  aa  (=  yaxa  );  activity  of  the  species  a  ,  here  ya  is  the 

activity  coefficient  of  the  species  a  to  account  for  the  nonideal  thermodynamic  behavior 
of  the  ‘membrane  solution  mixture’  ;  p,  the  thermodynamic  pressure  ;  ®  ,  the 
electrostatic  potential  [volts];  za  =  number  of  positive  or  negative  charges  on  a  particle 

(i.e.  number  of  charges  on  an  ion),  note  that  za  for  a  charge  neutral  species  is  zero;  R, 

the  gas  constant  =  8.314  J/mol-K;  T,  the  absolute  temperature  [K]  ;  F  =  Faraday’s 
constant  =  96487  Coulomb  per  gram-  equivalent,  i.e.,  96487  Coulombs  of  charge  on  1 
gm-mole  of  unity  charged  particles  such  as  H+  ions  or  e  electrons. 


Equations  (67)  for  the  various  mobile  species  are  rewritten  as: 

v 


allspecies  v  v 

Z 


p= 1 

*a 


'ap 


V  In  <7  +-^Wp  +  zo 

tl  rj-i  1  O 


f  F  ^ 

V® 

[rtj 

(68) 


(applied  for  a  =  H  + ,  CH3OH,  H2  O,  etc.) 

Note  that  terms  within  the  brackets  are  the  driving  potentials  for  the  transport  of  a 
species  a  through  the  membrane.  For  the  assumption  of  the  incompressible  solute  (H  + , 
CH  3  OH,  H2  O)  in  the  membrane  phase  and  system  (i.e.  membrane  phase)  isothermality, 
Eq.  (68)  can  be  written  as: 

zFO^ 


allspecies  v  v  / 

z  ^ 

P= 1  Vap 

*a 


V  -v 


-  X. 


lna„  +^-  +  — — 
“  RT  RT 


(69) 


(a  =H  +  ,CH30H,H20) 

The  expression  within  the  parentheses  plays  the  role  of  “total  potential”  or  driving  force 
for  the  transport  of  a  species  a  at  any  spatial  point  in  the  membrane  phase. 


The  general  boundary  conditions  are  set  up  by  requiring  the  “total  potential”  for  each 
species  a  to  be  continuous  across  the  boundary: 

At  x  =  0,  the  anode-side  boundary  of  the  membrane  in  contact  with  the  anode: 
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pva  zaF<£> 

In  a  H - 1 - 

RT  RT  @x=o,a 

{anode  side ) 


This  equation  is  simplified  to: 
u  .  \  (  v „  i 


pva  zaFO 
Inn  H - 1 - 

RT  RT  @x=o,m 

{membrane  side ) 

(for  species  a ) 


(]naa,oM  ~]na«,o,A)+  ~±  (Po,m~Po,a)+  -tp-  (®o^-®m)  =  0 


|  aa,o,M 
A 


+  [  ]  (po,M  -Pom)  +  ["^r]  K,m  -°m)  _0 

\K1  \K1 


(70a) 


(70b) 


v  ua,o,A  )  J  V  ^  J 

{a  =  H  +  ,CH30H,H20) 

A  simplification  results  if  the  assumption  of  chemical  equilibrium  across  the  anode- 
membrane  interface  is  valid. 


l.e.,  Oa  o  M  aa,o,A  • 


Po,a)+  l  (®o,M  0 


(a  =H  +  ,CH30H,H20) 


Also, 


v  |  z  F  (  v  i  z  F 

n  +— 2 _ (T)  =  _2_  n  +— " _ ® 

-t  o,M  fp  o, AA  rj-i  j  o,A  o,A 


(a  =H  +  ,CH30H,H20) 

For  the  one  (x-)  dimensional  case,  Eq.  (69)  becomes 
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(70c) 


Equations  (70b)  and  (70c)  may  be  written  as: 

ln«a,0,M  +{^\PoM  +^FF®o,M  =  ^aa0'A+[^^Po'A+^-<S>0,A  (70d) 

V  K1  /  K1  V  /  K1 


(70e) 


(71a) 


If  c  =  total  mixture  concentration  in  the  membrane  phase,  Eq.  (71a)  is  rewritten  as 


35 


allies  [CXa  )  U  ) 

2.  \Va-Vp) 


P=\  C~  ^aP 


d_ 

dx 


In  a  +  Fjl  + 


PVa  ,  ZaF® 


V 


RT  RT 


(71b) 
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(72) 


(a  =H\CH30H,H20) 

where  Na  =  molar  flux  of  species  a  along  x-direction  at  any  point  in  the  membrane 
phase. 

Now,  apply  Eq.  (72)  to  H  + ,  CH  3  OH,  H  2  O  : 

First  for  a  =  H+  : 
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Simplify: 
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Now,  first  we  write  the  flux  equations  for  methanol  and  water;  we  then  come  back  to 
Eq.(74a).  Say,  a  =  CH3  OH  in  Eq.  (72): 
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Divide  by  c 
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Now,  for  a  =H20  in  Eq.  (72): 
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At  this  stage,  it  is  important  to  note  that 


®  Pa 


(equivalent  to  the  Onsagar  reciprocal  equation) 


(75a) 
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4.1.3  Transport  Set  III  -  Flux  Equations  using  Onsager 
Thermodynamics 

4.1. 3.1  Introduction 

The  approach  is  based  on  the  multicomponent  diffusion  using  the  Onsagar 
thennodyanamics  as  given  in,  “J.S.  Newman,  Electrochemical  Systems,  p.265,  Prentice- 
Hall,  1991  edition.”  The  membrane  phase  is  assumed  to  be  “homogenous”  phase  with 
the  membrane  matrix  playing  the  role  of  a  “true  solvent”  for  the  permeating  species;  H+, 
CH3OH  and  H20. 


4.1. 3.2  Mathematical  Derivation 


The  starting,  multicomponent  diffusion  equation  for  the  transport  of  any  species  i  is 


Aft  =2X(v, -V,)=  RTZ-LJ-fa -v,) 


j  CT^ij 


(76) 


where  ci  is  the  molar  concentration  of  species  i,  cT  is  the  total  concentration  of  the 
species  including  the  solvent,  R  the  universal  gas  constant,  T  the  absolute  temperature, 
//,  the  electrochemical  potential  of  species  i,  K , .  the  friction  or  interaction  coefficients, 


DtJ  the  diffusion  coefficient  describing  the  interaction  of  species  i  and  j.  v,  is  the 

average  velocity  of  the  species  i  but  not  the  velocity  of  individual  molecules.  Equation 
(1)  is  similar  to  the  Stefan-Maxwell  equation  [R.B.Bird,  W.E. Stewart,  and  E.N. 

Lightfoot,  Transport  Phenomena,  p.570,  John  Wiley  and  sons,  Inc.  New  York,  (I960)] 
and  is  equivalent  to  that  developed  by  Onsagar  [Lars  Onsagar,  “Theories  and  Problems  of 
Liquid  diffusion,  “Annals  of  the  New  York  Academy  of  Sciences,  46  (1945),  pp.241- 
265]. 

The  molar  flux  of  species  i  is  given  by 

=clvl  (77) 


The  drag  coefficients  K , ,  and  diffusion  coefficients  D, ,  are  related  as  follows: 

1 J  1 J 


K, 


RT  Cj  Cj 

CT  Dij 


Equation  (1)  can  also  be  written  as 


(78) 
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(79) 


-?/) 

j 

In  this  equation  Ki  j  (if  -  v . )  is  the  drag  force  exerted  by  species  j  on  the  species  i 

because  of  the  relative  velocity  difference.  The  sum  on  the  right-side  of  Eq.  (79)  is  the 
total  drag  force  on  the  motion  of  species  i  per  unit  volume.  The  left-side  of  Eq.  (79), 

cj  (-  V//, )  is  regarded  as  a  driving  force  per  unit  volume  acting  on  species  i.  This  force  is 

considered  to  cause  species  i  to  move  with  respect  to  the  surrounding  fluid.  Equation 
(79)  expresses  the  balance  between  the  driving  force  and  the  total  drag  exerted  by  the 
other  species.  By  Newton’s  third  law  of  motion,  i.e.,  action  is  equal  to  reaction, 


D( .  =  D  u  (Onsagar  reciprocal  relation) 


(80) 


It  is  noted  here  that  the  number  of  independent  equations  of  the  type  (76)  is  one  less 
than  the  number  of  species  and  adding  them  leads  to: 


Xt'X/z]  =  X  2X(vy-v,) 


(81) 


The  left  side  of  Eq.  (81)  is  zero  as  evidenced  by  the  application  of  the  Gibbs-Duhem 
equation  at  constant  temperature  and  pressure  [J.W.Tester,  M.Modell,  Thermodynamics 
and  Its  Applications,  p.  148,  Prentice  Hall,  (1997  edition)].  The  right  side  is  zero  because 


Equation  (76)  is  quite  general  because  here  the  driving  force  on  a  species  i  is  related  to  a 
linear  combination  of  resistances  opposing  its  motion  instead  of  just  one  resistance 
relative  to  the  solvent.  The  spatial  gradient  of  the  electrochemical  potential  acts  as  the 
driving  force  for  diffusion  and  electric  migration.  The  number  of  transport  properties 
Dij  defined  by  Eq.  (76)  is  \n(n  - 1)/ 2]  for  n  component  system,  because  Dj  j  =  Dji  and 

Du  is  not  defined;  where  n  is  the  number  of  species  involved. 

For  the  membrane  phase  involving  membrane  matrix  m,  H+,  CH3OH  and  H20;  there  are: 
4(4- 1)/2  =  6  transport  properties  of  type  Z)  . . 
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Following  Eq.(76) 

c?Mi=  ZM,/(F,  -vm) 

j 


(82) 


Where  vm  =  velocity  of  the  membrane  material  which  is  zero  in  the  current 
situation. 


X M-J  Vj=  C,V  fi, 

j 

i  =  H+,  M,  W;  j  =  H+,  M,  W,  m; 
Mij=Kij  ;i*j 

Mu  =Ku -  HKik  -,i  =  j-,k  =  U+,  M,  W,  m. 

k 

Applied  to  i  =  H+;  M  (  CH3OH);  W  (H20): 

TMH+jVj=CH^VH+ 
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TjMMj  Vj  =  CMV/JM 
j 

YjMivj  Vj=cwV/iw 

j 

M  +  +v  ++M  +  v..+M  +  vw=c  +V  u  + 

H+H+  H+  H+M  M  H  W  W  H+  H+ 

M MH+  VH+  +ktMM  Etf  +  M Mw  Vffr  +  M Mm^rn  =CM^Mm 

M wh+V H+  +  Mwm  E M  +  kiWW  V W  +  MWmVm  =cw^MlV 


(82a) 


(82b) 

(82c) 

(82d) 

(83a) 

(83b) 

(83c) 


Note:  No  need  to  write  an  Eq.  for  membrane  (m)  species,  because  only, 

(n-1)  =  4-1=3  independent  Eqs.  of  the  type  (76)  exist.  Here  vm  (membrane  velocity) 
is  zero  (since  membrane  is  assumed  to  be  stationary  with  regard  to  the  laboratory 
reference  frame  at  the  steady  state  conditions).  So,  Eqs. (83a)  to  (83c)  become 

Mh+h+  V  +Mh+m  Vm+Mh+w  vw  =ch+Vmh+  (83d) 

M  mh+  ^ h*  +  ^mm  v m  vw  =cJmm  (83e) 

M WH+  V H*  +  ^WM  Vm  +Mww  Vw  =  cwVfiw  (83f) 
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K«  *K„;  *K»J 


2X, 


=  K, 


\K  a- K  A- K  A- K 

I^MH+  '  ^  M  M  '  ^  M  W  '  ^  M  m 


Mmm~  m  h>  +  ^  MW  K  M  m  ) 
^H*  M  +  K M  W  +  Km  m  j 

i=W(H20 ): 

rix 


Mww  Kww 


"  W  k 


Kww  WH *  WW  WM  ^Wm  ) 

M WW  =  ~{^WH*  +  ^WM  K-Wm  ) 

H*W  +  ^WM  +  ^Vm  ) 

Note:  Total  interaction  coefficients  involved  are: 


(84a) 

(84b) 

(84c) 

(85) 


(85a) 

(85b) 


(85c) 
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K H +  M  ’  K H +  ^  ,  Km  w  , 

^  II '  m  ,  .V/  m  ,  (On  * 


their  #  is  =  n(n-l)/2  =  4(4- 1  )/2  =  6 


Three  Eqs.  (83),  (83e),  and  (83f)  are  rewritten  as 


M  +  +  M  +  M  + 

//+//+  //+M  //+(f 

M mh*  ^  MM  M m  W 

M fVH+  M  WM  Mww 


Ch^Mh 

Cm^Mm 

cwV/Uw 


Call  the  matrix  made  of  M’s  as  L  matrix,  vM  column  vector  as  V  ;  and  the  one  in 


right-side  as  D  driving  force  vector;  Therefore  we  can  write  as 

L  V  =  D  (88) 

If  L  is  non-singular  matrix,  i.e.,  determinant  (L)^  0;  so  that  L  '(inverse  of  matrix  L) 
exists,  then  Eq.  (87)  or  (88)  is  solved  by  writing 

V  =  LXD  (89) 

(  Advanced  Mathematics  by  M.R.  Spiegel,  p.347  (1971  edition),  24th  printing  (1996)) 
The  system  of  equations  has  unique  solution. 

For  one  (x-)  dimensional  case,  Eq.(89)  is  written  as 


d  v,  =  /.  d  /) 


(89a) 


v  =  L  1  D 


(89b) 


where  vx  =  vu^x  ;  Dx  =  cM 


H  dx 
dMM 


(89c) 
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Note:  Find  L  ,  then  get  expressions  for  v H,  ^ ,  vM  x  and  vW  x  .  Then,  move  on 


Matrix  (L) 


L 


f  M  M  M  N 

M  MHi  M  mm  M  m  W 
M WH*  M  WM  Mww 


(90) 


Because  Mh+m=  Mmh+;  Mwh+=Mh+w  ;  M 

Note  that  the  L  matrix  is  the  symmetric  matrix. 
The  inverse  matrix  L  1  is  defined  as 

KJ 


=  M 

WM  M  W 


( L)  = 


det(z) 


(91) 


where  det  (L)  =  value  of  the  determinant  of  the  square  matrix  L,  and 
(M  jk)  is  the  matrix  made  of  cofactors  of  the  elements  in  the  matrix  L;  [m  jk  J  is  the 
transpose  of  [m jk). 

K)= 

(-D 
(-D 
(-D 


1+1 


2+1 


3+1 


^  MM  M  MW 
Mmw  Mww 


(-1)“ 


M H*  M  Mmw 
M  a*  w  Mww 


Mh+m 

Mh+w 

(- i)2*2 

Mh+h+ 

M  H+W 

MMw 

M 

1V±  WW 

M  H+W 

M 

1V±  WW 

Mh+m 

M  H+W 

(-D“ 

M H*  H* 

M  H+W 

M  mm 

Mmw 

Mh+m 

M  MW 

(-1) 

(-1) 

(-1) 


2+3 


3+3 


M H+M  ^MM 

M  H+  w  M  M  w 

Mh+h+  Mh+m 

M  H+w  M  MW 
M h*  H+  M h+ M 
M H*  M  Mmm 


Kd= 

(m  M  -M  M  )  -\M  M  -M  M  I  \M  M  -M  M  I 

V  M  M1V1  WW  1V1MW1VJMW)  V  H+ M  WW  H+ W  MW)  V  H+  M  M  W  H+  W  M  M ) 


'(» 


-\M  +  Mww  -M,  1WM  ,+ 

H  M  WW  MW  H  W 


)  (m  +  +MWW  -M  +  M  +  )  -(m  +  +M -M  +  M  +  ) 

/  V  //+//+  V  H  H  MW  H  M  H  w ) 


(m„-mmmw -mh,wmuu)  ~(m . Mww -M . U, 


H+  H+  M  W  H+M  H+  W 


)  (m  +  +Mmm  -M  +,  M 

J  \  U+  U+  M  M  U+  \4  l 


H+  H*  MM  +  '+  h+M^  H+M  JJ 
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K)= 

(m  M  -M2  )  \M  M  -M  M  I  \M  M  -M  M  M 

K  M  M1V1  WW  1V1  MW  J  \V1H+W1V1MW  H+ M  WW  J  V  H+  M  M  W  H+  W  M  M  J 

-<•  J, 

KJ  = 

(m  M  -M2  )  \MM  -M  m  )  (m  M  -M  m  )V 

V  M  M1V1  WW  1V1  MW  J  V  H+  W  MW  1V1  H+  M  WW  ]  \  //+  m  //+  W  MM  ] 

K-)r  = 

(m  M  -M2  )  \M  M  -M  M  I  \M  M  -M  M  I  ^ 

XV1  M  M1V1  WW  1V1  MW  J  \JKi  MW1V1  H+  W  1V1  H+  M  WW  )  \ LV1  H+  M1V1  MW  1V1  H+  W1V1  MM) 

-M„.„Mwr)  (m, -Ml„) 

I 

(92) 

(M„  Mn 

(M,Jr=  M2,  M22  M22  (93) 

vA/31  4f32  Af33  y 

where  Mn  —  (mmmMww  —M~mw^, 

M\2  —  \MmwMh+  w  —  M H+  mMww  j, 

M\i  —  [m N+ mM M  w  ~ M WM M M  j, 

M-  2\  —  H+  WM uw  —MH+  mM ww  j, 
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^23  h+  mM  H+w  M  H+ h+M  mw), 

M$\  —  [m h+  uM mw  —  M H+  WM M M  j  , 
Mn  =  (Mh+mMh+w-Mh+h+Mmw), 
=  (M  +  +Mmm  ~M\  ). 

33  V  H  H  MM  h*  M  ) 


The  determinant  of  matrix  (L)  is 

I M 

Det  (L)  = 


H*  H* 


M M  h+w 


mh +  Mmm  Mmw 


M. 


W  H* 


The  determinant  det  (L)  is  the  3rd  order  determinant . 


M wm  Mww 


From  Eq.  (91),  (93),  the  inverse  matrix  ( L  1 )  is 


'Mn 

Ma 

M„' 

M21 

M  22 

Mn 

Mi2 

j 

det(z) 


1 

(Mn  Mn 

4T21  4/22 

4T32 

Mu^ 

M2i 

j 

det(z) 

Mn 

Mn 

Mu 

det(z) 

det(z) 

det(z) 

M21 

4^22 

m22 

det(z) 

det(z) 

det(z) 

M31 

M32 

det(z) 

det(z) 

det(z) 

We  now  use  Eqs.(89b),  (89c)  and  (96): 


(94) 


(95) 


(96) 
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<  M„ 

m12 

(  dpi  ' 

fv  1 

det(z) 

det(z) 

det(z) 

<A+  , 

dx 

VM,x 

m21 

M  22 

M23 

r  dPiu 

det(z) 

det(z) 

det(z) 

CM  j 

dx 

Vw,x  ) 

M31 

M32 

M33 

dflyy 

v  det(z) 

det(z) 

det(z) ) 

r 

y  dx  J 

yM,x 

yVW,x  J 


(97) 


Hence,  equations  for  the  “net”  velocities  of  species  H+,  M  (methanol,  CH3OH),  and 
W  (water,  H20)  are  given  below: 


H* 


VM,X  = 


M, 


A 


det(z) 


dp. 


M 


21 


det(z) 


(f 


M 


A 

y 

c 

A 

V 


dx 


M 


12 


det(z) 


'M 


dfi, 


dx 


+ 


M 


22 


det(z) 


A 

y 

A 


dyM 

dx 


+ 


M 


13 


det(z) 


dfii 


"M 


31 


det(L) 


A 


d/u„  ,  \  M  v 


dx 


+ 


-  32 


det(L) 


A 


dx  j 

dMu 

dx 


M 


23 


V 

a  r 


det(z) 


+ 


M 


A 

y 

c 

A 

V 


djuw 

dx 

d fl  ]] 


7; 


33 


det(L) 


A 


dx 

dpi  if 

dx 


A 


J) 


Note:  det(L)  from  Eqs.(90),  [84a,  b,  c;  85a,  b,  c] 
M’s  from  Eq.(94),  [84a,  b,  c],  [85a,  b,  c] 
We  may  write  Eq.(98)  as: 


1 


H*,x 


det(z) 

1 

det(z) 


AV 


y  y  dx  j 


+  M 


12 


dpi  a 


dx 


M 


21 


dpi. 


v  H  dx  j 


+  M 


22 


dju 


M 


dx 


+  M, 


+  M 


dpp 


A 


dx 


A 


23 


d/u 


A 


w 


dx 


A 


(98) 


47 


"W, : 


dct(L) 


M 


31 


djUr 


v  H  dx  j 


+  M 


32 


dju 


M 


dx 


+  M 


33 


djuw 

dx 


Treatment  of  chemical  potentials  //,  of  the  species: 
For  a  charge  neutral  species, 

Vi=ni(@T,P)  +  RT\n(rixi) 


(99) 


(100) 


where  //,  =  chemical  potential  of  a  charge  neutral  species  i  in  the  membrane  phase  at  any 
point, 

Q:(@T,P)=  chemical  potential  of  species  i  as  pure  at  the  temperature  T  and  pressure  P 

and  in  the  same  state  of  aggregation  as  that  of  it  in  the  “membrane  solution 
phase”, (may  be  assumed  to  be  the  liquid  phase); 
yi  =  activity  coefficient  of  species  i;  it  is  a  function  of  temperature,  pressure  and 
composition  of  the  “membrane  solution  phase”; 
x,  =  mole  fraction  of  species  i  in  the  solution  phase. 

Note:  Here,  i  =  CH3OH,  H2O,  and  the  membrane  material  m  if  it  is  assumed  charge 
neutral. 

From  (100), 

djui  =dD.t .+  d[RT\n(yi  x,.)]  (101) 

=  (VtdP  -  StdT)  +  d[RT  ln(y,  xt )]  ( 1 02) 

dfj,,  =  V'dP  -  S^T  +  R  d[T  In {y.  x, )]  (103) 

where  Vt ,  S,  =  molar  volume  and  entropy  of  pure  species  i  at  (T,P)  of  the  system  at  a 

spatial  point  in  the  membrane  phase  solution. 

If  the  condition  of  isothermality  is  applied;  Eq.(103)  becomes 

d/u ,  =VtdP  +  Rd [T In (y.  x, )]  ( 1 04) 

divide  by  dx: 


'dMi\_v  dP  |  RTdln(y,.xi) 
v  dx  )  '  dx  dx 


(105) 


48 


i  =  M  (CH3OH),  W(HoO),  m  (membrane  material  as  a  component). 
Electrochemical  potential  of  H+  is  defined  as: 

f**=f*2m+zrr+FO 


RT  In 


lO.chem 

V  rt 


m , 


m 


+  zh.FO 


where  m°  =  a  reference  molality  (e.g.  1  g-mole  per  kg  of  a  solvent). 
The  molality  m  +  is  given  also  as: 


(106) 

(106a) 


mH:  =  H+  (106b) 

c0M0 

c0  Mo  =  solvent  concentration  in  mass  units,  here  solvent  assumed  to  be  mixture  of  the 
membrane  material  m,  water,  methanol;  for  the  transport  of  H+  ions.  M 0  is  the  molecular 
weight  of  the  solvent. 

The  expression  for  mH+  is  rewritten  as 


[chJct\ 


TYl  —  — — — 

M+  f  A 


\CT  J 


VXo  Mo  J 


(106c) 


where  cT  =  total  mixture  concentration  and  c’  =  hydrogen  ion  concentration. 


One  may  then  say,  x0  =  1-  x//+  ,  if  all  other  species;  membrane  material  (m),  CH3OH, 


H2O  are  taken  to  constitute  the  solvent  for  H+  ion  transport. 


Putting  Eq.  (106c)  into  (106): 


juH+  =RT  In 


lO.chem  „ , 

V  r,. 


H  H  xMm° 


+  zh+F  O 


[RT\nX\icliem\  +  RT  In 


rH+  — 

V  J 


RT  In  (M0m°)+  zH,F<& 


(106d) 

(106e) 
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dju. 


dx 


d 


R  —  [t  In  A0r’c+hem]  +  R 
dx  dx 


T  In 


7  H 


H+ 


V 


J 


0  ,  /, ,  o\dT  dO 

RlnlM/u  —  +  zij+F - 

'  dx  dx 

(106f) 


Under  the  condition  of  constant  temperature  and  assuming  that  the  effect  of  pressure 


variation  on  A,°'c+hem  across  the  membrane  for  the  pressure  variations  of  the  order  of  10  bar 


H 


is  negligibly  small;  then  Eq.(106)  reduces  to: 


d 


r djuu,  \ 

H  '  =  0  +  RT  — 
dx 


dx 


J 


In 


7  H 


H * 


V 


J 


0+zh+F^- 

dx 


d/u. 


dx 


dnH A 

v  H  dx  j 


RT  — 
dx 


=  E,.*7') 
=  Md 
=  E,.*r) 


- 

In 

v  + 

7  +  — — 

+  z  +F - 

H  dx 


d_ 

dx 


- 

f  „  \ 

In 

X  LJ  + 

7  + 

'  H+ 

l  )\ 

+  (w)F 


dO 

dx 


d 


dx 


In  y  +  +  ln- 


v 

d  In  y 


-  +  - 


d 


dx  dx 


J 


,  XH+ 

In  -J±- 


+  (w)F 


dO 

dx 


V  j 


dx 


From  Eq.(30) 


fcA0 

=  Cj  vt  —  +  c;i?r 

dilnYi)  |  d (lnx, ) 

v  dx  j 

dx 

dx  dx 

(106g) 

(106h) 


(106j) 


where  i  =  CH3OH,  H2O. 

Assuming  y u<  ,  7,-  variation  with  x  over  the  “membrane  phase”  negligibly  small, 
dlnyH+  _  d\nyi 


dx 


dx 


dx 


0;  then  Eqs.(106i)  and  (106j)  becomes: 

,  /  \  77  dO 

V  H  H  ’  dx 


c  +  CjRT——  [in  X;  ] 

llj  l  1  L  l  J 

dx  dx 


(  ^ 

=  (r 

"  d 

(  x  A 

In  H 

C 

v  H  dx  j 

\ch+K1  / 

dx 

111 

(107a) 

(107b) 
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Note: 


r  mu  1 

1 

m5 

det[z] 

J.S 

J.S 

in 


So  that 


r  mu  ~ 

(  dp  A 

m 5  J 

m 

det[z] 

V  ‘  dx  j 

J.s  m 4 

_  s  _ 

Simplification: 

dP 

If  — =  0,  if  pressure  on  both  sides  of  the  membrane  are  same  or  if  the  effect  of 
dx 

the  transport  of  a  species  is  negligible;  then  (107b)  becomes: 


dp, 

dx 


df1  ii 
dx 


c,RT 


=  cmRT 


d  lnxj 
dx 

i  =  CH3OH,  H20. 
,  d  lnx,. 


=  cmRT 


dx 

Jinx,,  dx. 


dxM  dx 


{xmct){RT).  1  d%M 


xM  dx 


(c ctRT )^- 
KT  ’  dx 


Also,  cu,  ^ Hl°  =cmRT  d  ln 


dx 


xwcrRT 


dx 

d  ln  x,„  dx,. 


Vyy  WJVyy 


dxw  dx 


=(xwctRT)  1  dx‘r 


xw  dx 


(< ctRT ) 


f  dxw  ^ 


v  dx  j 


Als°,  c  ^-[ln(xH+  /xa)]=  (xtftcr)y[ \n{xH+  /xj 


d 


dx 


dx 


=  (xH+Cr) 


d ln(x H+  / x  )  d(xHJxa) 
d[xH+  /x0)  dx 


dP 

—  on 
dx 

(108) 

(109a) 

(109b) 
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Substitute  from  (109a),  (109b),  and  (111a)  into  Eqs.(98): 


The  molar  fluxes  of  the  species  are  given  as: 


(112a) 

(112b) 

(112c) 

(113) 
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The  continuity  equation  for  a  species  i  in  the  membrane  phase,  under  the  steady-state 
condition, 

dN. 

— ^  =  0  (114) 

dx 

(i  =  H+,  CH30H  (M),  H20  (W)) 

leads  to: 

NiiX  =  Nix=0  =  Nix=Sm  =  const.  (114a) 

Nj  x=0  ,  Nj  x=s  are  the  fluxes  at  the  end  of  the  membrane  but  on  the  “membrane  phase 


side”. 


(116b) 
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N, 


w„ 


(Ur,*cir) 


CTXWVW,x  XWCT 


'  «»  v 


dct(L)  j 
M 


/  \  dX‘  jj+  dO 

{ctRT)——  +  zh+ctFx  h+  — 
y  dx  dx  j 


+ 


-  32 


det(z) 


(< ctRT ) 


dx„  W  M  v 


A 


dx 


-  33 


det(z) 


A 


(ctRT)A-' 
dx 


N, 


W  ,x 


M3\ 

det(z) 


(c2TRT) 


dx. 


dx 


+ 


M 


32 


det(z) 


[c2tRt) 


dx 


M 


VW 


dx 


+ 


(117a) 


M 


33 


det(z) 


[c2tRt) 


dx 


xn 


w 


dx 


+ 


M 


31 


det(z) 


(117b) 


y  uv  l  yx-/  y  y  y  wa  y  y  v 

Using  Poisson  equation,  it  was  found  for  the  membrane  phase  in  ‘Transport  Set  of 
Equations-I’: 


V  dx  j 


-  o 4  -  o c 


+ 


VC~**FV„-  2x) 


2  e 


(118) 


J 


where;  c  +  =  average  concentration  of  H+ with  the  entire  membrane  phase  taken  as  the 

mixture  of  [mobile  H+  in  the  fluid  phase],  CH3OH,  fUO  and  membrane 
material. 

=  XhXt 


8m  =  thickness  of  the  membrane  phase. 


F  <70 

© 

1 

e 

0 

!■  + 

r  zh+xh+Ctf ) 

y  dX  J 

Sm 

l  J 

(Sm  -2x) 


Putting  the  above  equation  into  Eqs.  (115b),  (116b),  (117b): 

N„+  = 


(119) 


f  Mu  A 

det(z) 


(c2RT) 


1  dx2H+ 

2  dx 


+ 


Mu 

det(E), 


'  Mu  V 

,det(z)A 

{zh+F  ct) 


( c2rRT  )xf 


dxM  H  M  v 


dx 


-  13 


det(z) 


x 


H * 


-\oA  -o 


c\  r 

—  + 


A 


z„+x„+ctF 


£ 


2  e 


"  K-  2x)\ 

(120a) 
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(120b) 


(120c) 


The  flux  equations  (120a),  (120b),  and  (120c)  are  quite  general  in  that  they  can  be 
applied  to  a  membrane  made  of  any  material,  provided  the  information  on  the  values  of 
the  interaction  coefficients  M ij  is  available.  The  following  equation  for  the  voltage 

distribution  in  the  membrane  is  from  the  solution  of  the  Poisson  equation  for  the  “net 
charge”  due  to  H+  ions  passing  through  the  membrane  in  the  company  of  H2O  molecules 
as  well  as  CH3OH  if  CH3OH  molecules  are  also  permeating. 


The  expressions  in  Eqs.  (120a),  (120b),  and  (120c)  indicate  the  fluxes  of  the  species  H+, 
CH3OH  (M)  and  H2O  (W)  via  diffusion  mechanism,  where  the  diffusion  flux  of  each 
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species  is  affected  by  the  gradients  of  the  concentrations  of  all  mobile  species,  and  via 
electric  field  effect. 


Equations  (120a)  through  (120c)  are  now  transformed  as  follows: 
For  the  treatment  of  Eq.(120a): 


dx 


dx 


d  ( 
—  \x 


dx 


d  i 


XM 


Xw 


dxH, 

dx 

dxH+ 

dx 


Using  the  above  equations,  Eq.(120a)  is  rewritten  as: 

N.,+  = 


r  m  ^ 

ivi  n 


det(z) 


(  c2tRT^ 


dx ; 


dx 


+ 


J  J 


det(z) 


( xH+xM )  xM 


dx. 


dx 


dx 


+ 


My 

det(z) 


M  n 
det  {L)y 


{ctRt\  ( xH+xw )  xM 


dx. 


dx 


dx 


-  QA  - cpc 


Of 


+ 


2  e 


J 


+ 


For  the  treatment  of  Eq.(120b): 


dx 


H * 


dx 


^-(xMxHl)-x 


dx 


M 


dx 


H*  /  H* 


dx 


dxu 


VM 


dx 


d  /  \  _  dx  yi 

——  \XM  xw )  —  xw  - 
dx  dx 


Using  the  above  equations,  Eq.(120b)  is  transformed  into: 


N 


M,x 


M 


21 


{c}RT)  y(%v)-X 


dclff-l 


dx 


dx 

ua,m 

dx 


+ 


M 


22 


det(z) 


f  c~RT^ 


J  L 


dx 


M 


dx 


M 


23 


(ctRt)  (xMxr)  xw 


det(l)y 


dx. 


dx 


dx 


M 


21 


det(E)y 


(vF4)(x«xfft) 


-0-0 

Of. 


+ 


■  + 


z  +c  +FY  . 

"  "  (S„  -  2x) 


2  e 


J 


(121a) 


(121b) 
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For  the  treatment  of  Eq.(120c): 
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Putting  this  information  into  Eq.(120c): 
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(121c) 


Equations  (121a),  (121b),  and  (121c)  describe  the  fluxes  of  the  species  H+,  CH3OH  (M) 
and  H2O  (W)  through  the  membrane  via  diffusion  and  electric  field  effect.  Note  the 
transport  of  CH3OH  (M)  and  FEO  (W),  via  electric  field  effect  also,  for  these  charge 
neutral  species. 

These  equations  are  quite  general.  Also,  for  the  condition  of  steady-state  the  fluxes 
N H,  x  ,  NM  x  and  NWx  are  invariant  with  respect  to  the  distance  x  through  the 

membrane  phase. 
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In  what  follows  an  approximate  method  (i.e.  engineering  practice  method)  is  followed  to 
relate  the  fluxes  to  the  species  mole  fractions.  Here,  diffusion  interaction  related 


parameters 


f  M  A 


s  are  assumed  to  be  constant  over  the  membrane  phase  thickness, 


vdet(X)y 

8 m  .  Total  membrane  phase  molar  concentration  cT  is  also  taken  constant.  In  fact,  one 
may  take  the  average  values  of  these  parameters  at  the  membrane  ends. 


Treatment  of  Eq.(121a): 

Multiply  both  sides  of  this  Eq.  By  dx  and  integrate: 
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Now,  some  approximations  have  to  be  made 
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Equation  (122f)  becomes 


59 


s,„  _  C$)A  _(J)C  ( z  J7C  x 

f  2  ,  ZH*rLTXH* 


1R 


5 - U  ^pl(«.-2,)kfc- 


fe.)  -b'' -$c|+[ z",'FCr't"' 


=  feOH0"*  _a)ClJ  (122i) 


Equation  (122b),  now,  becomes 
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Check:  Units  : 


Sm\ 

Units  of  each  of  the  terms  :  (1),  (2),  (3): 

nr  f  mol  ^  1  _  f  mol 

J.s\  \  m  J  \_m\  \m  -  s 

Term  (4)  : 

m 5 
J.s 

Equation  (123)  is  dimensionally  correct. 

Treatment  of  Eq.(  12 1  b): 

Multiply  this  Eq.  by  dx  and  integrate  over  the  membrane: 


coulomb\  f  mol  V  f  volt 


J.s  -nr 


Divide  both  sides  by  8m  : 


(124) 
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NM  = 
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(125) 


Treatment  of  Eq.(121c): 

Multiply  (121c)  by  dx  and  integrate  tenn  by  term  over  the  membrane  for  x  =  0  — »x 
8,-. 


Nw,x  8m  = 
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Comment: 

In  Eqs.  (123),  (125),  and  (127),  on  the  right-hand  side  of  each  equations,  the  terms 
contribute  to  the  flux  of  the  concerned  species  via  the  multi  component  diffusion  and  the 
last  term  on  the  right-hand  side  represents  the  flux  due  to  electric  field.  In  the  case  of  M 
(CH3OH)  and  W  (H2O),  the  last  term  is  defined  here  as  the  electro-osmotic  flux. 


Equations  (123),  (125)  and  (127)  are  now  written  in  terms  of  concentrations  rather  than  in 
terms  of  mole  fractions: 
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Now,  the  species  H+,  CH3OH  (M),  and  H2O  (W)  concentrations  in  the  membrane  phase 
at  its  boundaries  with  the  electrodes  are  expressed  in  terms  of  the  species  compositions  in 
the  porous  electrodes  using  the  concept  of  phase  equilibria  between  the  membrane  and 
electrodes  with  regard  to  each  species. 


For  the  liquid  mixture  in  the  porous  anode  in  equilibrium  with  the  membrane  phase  at  x  = 
0,  with  regard  to  a  species  i,  the  condition  for  equilibrium  is 
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(131) 


A—e  m-eo 
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(i  =  CH3OH,  H20) 
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^  m-eo  ^ 
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(i  =  CH3OH,  H20) 


(132) 

(133) 


where  ju  °L  =  chemical  potential  of  the  pure  liquid  species  i  at  the  temperature  and 
pressure  of  the  mixture  in  the  anode, 

ju°m  =  chemical  potential  of  the  pure  liquid  species  i  at  the  temperature  and 
pressure  of  the  mixture  in  the  membrane  phase, 
xf~A~e  =  mole  fraction  of  species  i  in  the  liquid  phase  mixture  in  the  porous  anode 
end  in  contact  with  the  membrane, 

c"'  eo  =  mole  concentration  of  species  i  [mol  per  unit  volume]  in  the  membrane 
phase  in  contact  with  the  porous  anode  end  at  x  =  0, 

^  ^  M-eo  ^ 

c  =  1  mol  per  unit  volume;  used  here  as  adjusting  parameter  to  make  — - 

l  Co  V 

a  dimensionless  quantity, 

yf~A  =  activity  coefficient  of  species  i  in  the  liquid  phase  mixture  in  the  porous 
anode  to  account  for  the  nonideal  behavior  of  species  i  in  the  mixture, 
yj'"'  A  =  activity  coefficient  of  species  i  in  the  membrane  phase  to 
account  for  the  nonideal  behavior. 

The  standard  state  chemical  potentials  for  the  liquid  phase  in  the  anode,  /u"L ,  and  for  the 


membrane  phase,  ju°'n ,  are  assumed  to  be  identical,  i.e.  ju°L  =  , 

Eqs.  (131),  (132), and  (133)  lead  to 
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This  leads  to: 
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(135) 
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(i  =  CH3OH,  H20) 


Here,  s{ 


L-A- 


7, 


L-A 


c,m-A 


7i 


,  has  been  defined  as  the  liquid  solubility  coefficient  of 


component  i  in  the  membrane  phase  corresponding  to  the  liquid  phase  composition  in  the 
anode.  Note  that  c'"  eo  is  the  molar  concentration  of  the  species  i  in  the  membrane  phase 
at  x  =  0. 


Following  the  similar  procedure,  it  leads  to 

m—e  L—C-m  L—C—e 
C,  =  Si  X ■ 


(137) 


(i  =  CH3OH,  H20) 

where  c,"!~e  =  molar  concentration  of  species  i  in  the  membrane  phase  at  x  =  Sm  , 


sL-c~m 


7 i 


L-C 


c,  m—C 


7 , 


the  liquid  solubility  coefficient  of  component  i  in  the  membrane 


phase  corresponding  to  the  mixture  composition  in  the  porous  cathode  end 
in  contact  with  the  membrane  at  x  =  8m  . 


If  the  dependence  of  the  solubility  coefficient  of  a  species  i  in  the  membrane  phase  on 
compositon  can  be  relaxed;  then,  one  may  set 

sf~A~m  =  sf~c~m  =  s'"  (solubility  coefficient  of  the  species  i  in  the  membrane 

phase)  (138) 


Equations  (136)  and  (137)  are  rewritten  as 
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where  s"‘  =  solubility  coefficient  of  the  species  i  in  the  membrane  phase  at  the  prevailing 


temperature  and  pressure  conditions. 


Treatment  of  H  : 


, .  A-e  _  m,eo 
" H+  " H+ 


(140) 


where  /uY+e  =  electrochemical  potential  of  H+  in  the  porous  anode  end  in  contact  with  the 
membrane  phase, 

=  electrochemical  potential  of  H+  in  the  membrane  phase  end  in  contact 
with  the  anode  at  x  =  0. 


The  electrochemical  potential  of  H+  is  expressed  as 


RT  ln(/L' 


(141) 


where  m/r  =  molality  of  H+  in  the  solution  in  a  porous  electrode  or  in  the  membrane 
phase, 

y  +  =  y H,  (temperature,  pressure  and  compositon)  =  activity  coefficient  of  H+, 

A? .  =  a  proportionality  constant,  independent  of  composition  and  electrical  state, 

but  characteristic  of  H+  and  the  solvent  and  dependent  on  temperature  and 
pressure. 


Equation  (141)  is  used  in  Eq.(140)  to  obtain 


nO,A  ,„„A-e  A-e _  nO,m  ^m-eo  „,m,eo 
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(142) 


If  it  can  be  allowed  that  =  A^+  (it  may  become  possible  if  one  can  assume  that 

water  is  solvent  for  the  mobile  H+  ions  in  the  porous  anode  and  membrane);  then, 
Eq.(142)  becomes 
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The  solubility  coefficient  of  H+  ions  in  the  membrane  phase  in  the  presence  of  water  may 
be  defined  as 
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where  c"{, eo  =  molar  concentration  of  water  in  the  membrane  phase  end  in  contact  with 
the  porous  anode  at  x  =  0, 

cf(! e  =  molar  concentration  of  water  in  the  porous  anode  end  in  contact  with  the 
water, 

Yh+  ’  y"‘r°  =  activity  coefficients  of  H+  ion  in  the  porous  anode  solution  and  in  the 

membrane  phase  solution,  respectively,  at  the  membrane-anode  interface 
at  x  =  0. 


Note  that  sLH™  A  is  controlled  by  the  product  of 
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From  Eqs.  (143b)  and  (144), 
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Following  the  similar  procedure  it  can  be  shown  that 
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where  c"  ,e  =  molar  concentration  of  H+  in  the  membrane  phase  end  in  contact  with  the 
porous  cathode  at  x  =  Sm , 

ccH:e  =  molar  concentration  of  H+  ion  in  the  mixture  in  the  porous  cathode  end  in 
contact  with  the  membrane  at  x  =  Sm  , 
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membrane  phase  at  the  cathode-membrane  interface. 
If  it  can  be  allowed  that 

L-m-A  _  L-m-C  _  e  m  . 

S  H+  S  H+  ^H+  9 

then,  Eqs.  (145)  and  (146)  are  written  as 
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Molar  concentrations  of  H+,  CH3OH  (M),  and  H2O  (W)  appearing  in  the  membrane 
phase  in  Eqs.  (128),  (129),  and  (130)  are  expressed  in  terms  of  the  concentrations  of  the 
species  in  the  porous  anode  and  cathode  using  Eqs.  (139a),  (139b),  (147)  and  (148)  to 
obtain 
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5  Conclusions 


Derivation  of  the  transport  flux  equations  for  a  direct  methanol  fuel  cell  solid  polymer 
electrolyte  membrane  using  i)  moderately  dilute  solution  theory  based  on  Fickian 
diffusion,  ii)  concentrated  solution  theory  based  on  the  generalized  Stefan-Maxwell 
equations,  and  iii)  concentrated  solution  theory  based  on  the  Onsager  irreversible 
thennodynamics  approach  to  transport  processes  have  been  presented.  The  developed 
equations  account  for  the  effect  of  the  voltage,  pressure  and  the  species  concentration 
gradients  on  their  transport  fluxes. 

These  equations,  describing  species  transport  fluxes  through  a  solid  polymer  electrolyte 
membrane  are  to  be  coupled  with  the  equations,  yet  to  be  developed,  describing  mass 
transfer  through  the  electrode  porous  backing  layers  and  species  mass  transport  with 
electrochemical/chemical  kinetics  in  the  porous  electrode  layers.  The  coupled  equations, 
describing  the  various  phenomena  occurring  in  a  DMFC  in  operation,  will  be  of  immense 
significance  in  the  evaluation  of  performance  and  design/development  of  a  DMFC. 

The  developed  transport  equations,  in  this  report,  can  be  used  to  experimentally 
determine  the  transport  properties  of  a  membrane,  from  very  dilute  to  very  concentrated 
solutions,  such  as  permeability  of  the  fluid  mixture  and  effective  mass  diffusivities  of  the 
various  species. 
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